PPP bP PP PP PP PP PP PP 


Oo oo N HY FW N y 
. 


10. 


Bh we eee ce ee eee 1 
基本 方程 eee 2 
. 粘性 应 力 和 热 通 量 Ne 10 
» 基本 方程 的 积 从 形式 pe 15 
. 相似 性 和 流动 参数 pe 19 
理想 气体 24 
. 理想 气体 的 非 绝 热流 动 ， 环 量 与 涡 量 ee 26 
。 理想 气体 的 绝热 流动 。 MARI BE rrerrrrererereresrrerresrerrereres 34 
无 旋 流 。 ia AE none tt 上 36 
无 旋 流 的 变 分 潜 ee 39 
绝热 定常 流 ene cece tee eee ten etetene ce cee eee eee nent tteneenenseeeeteeeae een ees 41 
. 二 维 定 常 均 能 流 。 流 蜀 数 46 
. 二 维 定常 无 旋 流 。 速 庶 面 变换 51 
.对 完 双 气体 律 的 偏 窟 cee ceteeteee cee eenereeeensceeeeesseeseesetan eee se ees 55 
. 一般 正 交 坐标 季 中 的 运动 方程 的 表达 式 eet ece eter c ete rere eee eeees 58 
郊 考 文献 ee eee eee eer 66 


Al. 引 8 


对 设计 航空 器 或 具有 流动 工 质 的 机 械 的 工程 师 来 说 ， 他 所 需 
要 的 资料 是 运动 流体 作用 在 一 个 特定 物体 或 多 个 物体 表面 上 的 压 
力 、 切 应 力 、 温度 及 热 通 量 矢 ; 这 些 物体 是 处 于 在 特定 条 件 下 流动 
着 的 流体 之 中 的 。 气体 动力 学 的 主要 任务 就 在 于 提供 这 些 资料 . 
除 上 述 的 资料 外 ,再 加 上 流体 的 速度 矢 和 密度 这 两 项 资料 ,一 般 地 
认为 有 关 流 体 状 况 的 资料 就 齐全 了 。 所 以 在 指定 空间 域内 各 点 
上 ,在 指定 时 间 范 围 内 每 一 瞬间 ， 上述 各 量 一 旦 确定 的 话 , 那 末 气 
体 动力 学 问题 就 算是 完全 解决 了 .本 统 的 目的 在 于 系统 地 国明 一 
般 问 题 ,大 提供 作 分 析 所 需 的 诸 方 程 .- 

不 过 在 进行 仔细 分 析 以 前 ， 最 好 先 计算 一 个 在 气体 动力 学 中 
十 分 重要 的 量 , 那 就 是 小 扰动 的 传播 速度 ， 这 个 速度 常 称 为 声速 . 
由 于 在 流动 流体 中 的 物体 对 流体 产生 扰动 ,可 以 料 到 , 流 场 的 特性 
会 在 很 大 程度 上 决定 于 平均 运动 速度 与 声速 之 比 。 事 实 上 ， 仅 根 
据 这 个 比值 就 可 以 把 流动 问题 划分 为 四 类 : 一 、 亚 声速 流 , 郭 流体 
速度 小 于 声速 ; 二 、 跨 声速 流 , 朗 流 体 速 度 和 声速 差不多 大 小 ; 三 、 
超声 速 流 ,部 流体 速度 大 于 声速 ; 最 后 是 高 起 声速 流 , 郎 流体 速度 
大 大 超过 声速 . 

现在 来 计算 小 扰动 的 传播 速度 : 考虑 一 从 左 向 右 的 均 与 流 . 
这 流动 越过 一 微弱 突 跃 面 , 突 跃 面 位 填 不 动 , 于 是 流 场 是 定常 的 ， 
也 就 是 说 流 场 不 随时 间 变 化 。 流体 流 过 突 跃 面 时 ， 速 度 由 * JH 
u + du, 密度 由 2 HH po + do, 压强 由 Pp 增 为 p 十 4p. 现在 考虑 
一 具有 单位 截面 面积 的 流 管 。 在 越过 突 跃 面 时 ， 流 体质 量 不 能 有 
PIG IR, ik 

pu = (p + dp)(u + du), (1.1) 

又 越过 突 跃 面 时 ,动量 的 增 量 必 为 作用 在 流体 上 的 压强 增 量 所 斑 


ee HES: 


衡 , 故 
(p+ do)(u + du? — pu? = —dp, (1.2) 
车 仅 取 一 级 微量 , 则 从 方程 (1.1) 和 (1.2) 得 


显然 ,在 一 个 和 流体 一 起 运动 的 观察 者 看 来 , 流体 是 不 动 的 , 而 扰 
BASE u 传播 . 因此 , «确实 就 是 小 扰动 的 传播 速度 ,也 就 是 声 
速 。 合 声速 为 a, 则 

ad 

a 一 do’ (1.3) 
只 要 知道 压强 和 密度 之 间 的 关系 ， 就 可 以 从 方程 (1.3) 计算 出 声 
速 。 不 可 压 流 体 的 声速 是 无 限 大 ,这 是 因为 它 的 压强 可 以 变化 ， 
而 Pp 不 能 变 的 绿 故 ， 由 于 扰动 传播 不 能 产生 或 消灭 能 量 ， 所 以 这 
个 传播 过 程 必 是 绝热 的 ， 事实 上 ,一 种 特有 的 压强 -密度 关系 是 等 
箭 关系， 对 于 完全 气体 来 说 ， 邹 对 于 成 分 分 子 闻 的 相互 作用 可 以 
完全 忽略 的 气体 来 说 ,压强 对 密度 的 变化 率 是 

A =r, (1.4) 

式 中 的 7 是 比 热 比 ， 在 普通 温度 和 压强 之 下 ， 大 多 数 气体 是 完全 
气体 。 这 时 声速 的 计算 十 分 简单 , 朗 


a@ayl, (1.5) 
p 
A.2. 基本 方程 


AY — BAYER REE x, 此 处 i = 1,2,3, FEN RI 
t, 坐标 x; 点 的 流体 速度 分 量 记 为 wi, i 一 1, 2, 3。 其 密度 和 温度 
分 别 记 为 2 ST, 流体 的 温度 是 指 , 用 一 种 滞后 时 间 可 以 忽略 的 
天 和 流体 一 起 运动 的 温度 计 所 测 得 的 温度 ,全 qi(i = 1, 2, 3) 为 
热 通 量 矢 , 郎 单位 时 间 流 经 单位 面积 的 热流 量 . 

流体 中 的 应 力 可 以 表示 为 一 微 元 流体 各 个 表面 各 方向 受 力 的 
强度 。 图 A.2 a 示 边 长 为 dn, dm, dxs 的 立方 微 元 体 上 所 具有 的 
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O33 
33 十 —— dx; 
ax T32 下 gay dx) 
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DT23 
TW 
23 十 Er dx2 


Drr22 
m2 十 Ox dxz 
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图 A2a 
这 些 应 力 。 具体 说 来 ，ma 指 作用 在 与 n RM LMA 
' 轴 的 单位 面积 上 的 拉力 ,或 说 拉 应 力 。xw 是 作用 在 与 wi HE 
址 的 面 上 而 指向 为 = 轴 的 单位 面积 上 的 剪 切 力 , 或 说 剪 应 力 。 Br 
以 , xz 代表 作用 在 与 x; 轴 垂 站 的 面 上 而 指向 为 x; 轴 的 应 力 , 这 九 


。 3 ù 


个 应 力 的 组 合 简 称 为 应 力 张 量 ( 见 第 IV 4*, B 编 )， 当 然 微 元 立方 
体 有 两 面 垂直 于 同一 坐标 轴 . 例如 与 轴 垂 直 的 在 处 有 一 个 
面 , Æ n 十 ax 处 还 另 有 一 面 。 如 果 应 力 是 坐标 的 画 数 , 并 且 在 
左边 的 面 上 各 是 Ry Wir, M13 的 话 ， 那 末 在 右边 的 面 上 应 为 ma 十 
(Oru/Ox1) dx, my + (Omn/Ox) dzi, m3 + (Omn/On) dn, 微 元 
立方 体 的 其 他 儿 对 面 应 是 同样 情况 . 
现在 来 计算 诸 应 力 的 合力 分 量 . 在 方向 的 分 力 是 
[= + (=u + a dx,)| dx, dxs 


+ [ma + (xa + ora zj dx; dx, 


x2 


+ | -= 十 (xa + Ona dzs)| dx, dx; 
Oxs 


— (om + Oxn + ?za ) dx, dx, dx 


Ox Ox, x3 


= ( > San) dx, dx, dx3 = Ons dx, dx, dx, 
7=12,3 Ox; 


aj Xi 

最 末 一 式 采 用“ 取 和 约定 ”(Summation convention)， 也 就 是 说 凡 
在 一 项 里 下 标 重 履 ， 那 末 该 项 就 对 该 下 标 求 和 **。 所 以 各 应 力 的 

合力 在 i 向 的 分 力 是 
Orji 
xi 
合力 的 诸 分 力 必 须 与 流体 的 彻 体力 和 惯性 力 平衡 ， 显 然 一 立方 微 
元 流体 的 彻 体力 和 惯性 力 必须 正比 于 其 体积 dr da dx3; 否则 将 无 
法 和 应 力 的 合力 保持 平衡 .当然 ,事实 上 这 个 条 件 是 得 到 满足 的 ， 
因为 彻 体力 和 惯性 力 与 微 元 立方 体 的 质量 成 正比 ， 而 其 质量 正好 
Æ p dxri dxz dx, 在 建立 此 平衡 方程 之 前 ， 先 考虑 力矩 的 平衡 关 
* 指 “ 高 速 空 气动 力学 及 喷气 推 进 ” (High Speed Aerodynamics and Jet Pro- 

pulsion) 丛书 的 第 IV 移 ; 余 类 同 编者 注 . 

* 在 工程 或 物理 文章 中 常常 要 用 到 Darbi 之 类 的 表达 式 , 为 了 方便 , 现在 一 般 都 
采用 这 样 一 个 约定 ,即将 求 和 号 省 去 , 而 凡是 一 项 里 同一 下 标 出 现 两 次 时 ,就 次 


理解 为 它 是 有 求 和 号 的 。 Wath; = abı + asb + asbs; ajbkrck = aj(bicr + 
pacs + pacs) 一 一 译 者 注 . 


dx, dx dx;, 


-or 一 一 一 一 


条. 仍 取 微 元 立方 体 ， 计算 各 力 对 一 条 三 行 于 x3 轴 的 轴线 的 力 
2E, 彻 体力 和 惯性 力 的 力 距 是 与 微 元 体积 和 dx 或 dx, 的 乘积 成 
正比 , 是 一 个 四 阶 微量 。 图 A.2 a RMN Cn, s x) RMA 
行 于 凡 轴 的 一 根 轴线 的 力矩 平衡 关系 式 如 下 : 
= (ap dx, dx3) dx, — (xy dx; dx) dx, + 四 阶 微量 项 。 故 得 
my 一 Ry 
或 一 般 地 表示 为 
Mij 一 Tri。 (2.1) 
从 力矩 硅 衡 得 出 来 的 这 个 结论 总 是 成 立 的 ， 并 表 朋 应 力 张 量 只 有 
六 个 独立 的 分 量 。 这 样 的 张 量 称 为 对 称 张 量 , 
通常 热力 学 压强 ?是 由 p,p 和 7 了 满足 流体 的 状态 方程 来 定义 
的 : 
p = plT, p), (2.2) 
这 个 状态 方程 是 在 静止 条 件 下 定 出 来 的 。 所 以 “粘性 应 力 张 量 ” 
Tj 应 定义 为 
Ti = m; + ôP» (2.3) 
式 中 的 8 EEP REAR 8, 定义 为 
64 = l, i =j; 
=0, ixj. (2.4) 
AM i AGN, rz 一 xy, 而 方程 (2.1) 双 说 明 Ty = Tj 所 以 zi 也 
是 对 称 张 量 。 无 粘性 流体 是 ry 恒 等 于 需 的 流体 。 关于 区 分 流体 
ARIETA P (hydrostatic pressure) 和 粘性 张 量 zi; 的 基本 物理 意义 将 
FTH. 
取 一 个 边 长 为 dx1, dx, dxs 的 立方 体 来 看 。 单位 时 间 流 出 立 
ARPA BREE 
(Sees + .om 十 Sous) dx dx dx3, 


Ox Ox, x3 
IRDA AE PSE, A PRN ERS ETE) Eb, Ep 


Opu:! , Opu Ses) (2 ) 
ZEEL p eR y SEBS) dx, dx, da, = —| <E )\ dx, dx- d 
( Ox; Ox, arg) 00° Bt | nk O82 atas 


eile ined TEREE SEEE Te 


所 以 质量 守恒 方 程 或 连续 性 方程 是 
oo 十 Opu; 
Ot Ox; 
这 里 又 用 了 取 和 约定 。 也 就 是 说 
ou _ Opus 
Ox; £=1,2,3 Ox; ` 
根据 类 似 计 算 办 法 ， 得 单位 时 间 单 位 体积 在 i 方向 的 动量 增 
量 是 


一 0， (2.5) 


Opu; + Opu iu; 

Ot Ox; | 
JEZI ARW eh T OCI RA Se FE D 
i. 作用 力 有 两 个 来 源 。 设 X; 代表 单位 质量 所 受 的 来 自 外 源 的 
力 的 分 量 , 那 末 加 速 力 的 一 部 分 就 是 X;。 另 一 部 分 则 来 自 应 力 张 
量 rjj. 的 tA HA 


2 eH) BF, 
根据 张 量 rr va Beton en DL iR. 单位 时 间 单 位 体积 的 i 


向 动量 滔 培 量 必 等 于 作用 于 单位 体积 的 i 向 力 ， 表 明 这 个 关系 的 
动力 学 方程 便 是 


Sets 2 Comins) = pX; + ri, (2.6) 
Ot Ox; 
把 左 侧 各 项 展开 ， 守 用 上. 弱 作 方程 (2.5) ARO OLTES 
Ou, p ti a yg, 4 Or 
Oe + uj Bz, X: + PE (2.7) 
如 引信 粘性 应 力 张 量 Tas 则 方程 (2.6) 和 (2.7) 变 成 
pui 8 Coun) = — OP Or 
az + an (pu;u;) Ox, + pX; +e, (2.8) 
及 
Ou ,Ou a 1 OP 4 x 4 A or, (2.9) 
9 Ox; 2 Ox; p Ox,” 


方程 (2.6) 至 (2.9) 是 动力 学 方程 的 四 种 形式 . 


在 方程 (2.7) 及 (2.9) 中 , 左 侧 是 随 流 体 运动 计算 出 来 的 流体 微 


困 的 加 速度 ， 也 就 是 流体 的 实际 加 速度 ， 为 了 证 明 这 一 点 是 正确 

的 ,可 以 拿 任何 一 个 量 ft, xi) 来 看 ,这 是 时 刻 上 RADA x; 的 图 

数 . 假设 观察 者 和 流体 微 团 一 起 运动 ,到 了 +: 十 eg, ere 

间 点 的 坐标 为 x; + udt, 所 以 对 于 这 个 观察 者 说 来 , 物理 量 f 随 

时 间 的 变化 率 是 

. [fz + dt, x; + 4;6t) — f(z, x;)] _ (9. 3 
im o an E tu Z) Í. 

AMC BEDERAK E D/Dt, 则 

D 0 ð 


一 一 + Wi. 
Dt Ot Ox; 


如 。 代表 单位 质量 流体 的 内 能 ， 则 单位 时 间 内 单位 体积 此 内 
能 的 增 量 为 


(2.10) 


De 
Dt” 
此 能 量 增 量 来 自 三 方面 : 第 一 ,来 自 象 吸收 辐射 热 ,化 学 反应 及 燃 
烧 的 外 加 热 。 记 单 位 时 间 内 加 给 单位 质量 的 热能 为 0， 第 二 , 来 
自 热传导 , 其 热 通 量 矢 4 实际 是 减少 微 元 立方 体 中 的 能 量 ， 故 单 
位 时 间 单 位 体积 的 热能 增 量 是 

_ ôq 

Ox; 

第 三 ,来 自 应 力 张 量 xi 对 流体 所 做 的 功 ， 兹 对 此 作 仔 细 分 析 : 

先 考虑 xz 轴 方 向 的 一 维 流 动 的 简单 情况 . 取 一 根 单位 截面 的 
流 管 ( 见 图 A.2b). 在 £ AIRETA n Bl ay 十 dey 的 位 
E. 在 稍 后 的 时 列 1 十 dz, 原 在 x 点 的 流体 移 至 r + u de 点 ,而 
原 在 x, + dx, FRA BB x, + dxi + [u + (Ou;/Ox,)dx,] de 
A. 原来 这 段 流 体 的 体积 是 ar 1， 现 在 变 成 了 [1 十 
(Ou;/Ox,) dz ] dr1* 1， 所 以 体积 的 变化 是 《9xVaxm) dx, dt, BFK 
的 功 等 于 拉 应 力 xu RARE, Bp 


Ou, 
— dx, dt 
Tu Ox, 1 . 


六 一 B 


图 A.2b 
单位 体积 的 流体 在 单位 时 间 内 做 的 功 就 是 
TU Ba 
XL 
在 三 维 流 的 一 般 情 况 下 ,还 有 类 似 的 二 附加 项 : xw( Bw2/6x2)， 
xa Ous/ Oxs). 事实 上 Gu/9xi 就 是 流体 微 困 在 x 方向 的 拉 应 变 
OT BY [a] AE AES, Ou2/ Ox, 是 x2 HAEREERE, Ouz/Ox3 是 x3 
Paar AVE BSE ES, 所 以 单位 时 间 单位 体积 所 做 的 功 应 等 于 各 应 
力 和 对 应 的 应 变 变 化 率 的 乘积 之 和 ， 这 里 面包 括 切 应 力 及 切 应 变 
变化 率 ， 问 题 是 切 应 变 率 是 什么 ? 为 回答 这 个 问题 ， 先 观察 表达 
流体 微 团 在 运动 中 的 变形 率 张 量 6wi/6x;。 这 个 张 量 可 以 分 成 对 
称 部 分 和 反对 称 部 分 : 
Ox, Ox. Ox; 
Ou; _ | Ot, Ou Ou 


Ox; Ox, Ox, Oxs 


Ous Ous Oua 
Ox, Oxz Ox; 

Ox, 2 \ Ox, Ox, 2 \ Ox3 Oxy 
2 Ox, Ox, Ox, 2 Ox3 Ox 
(和 aa + 2u) 4( Sn + Ses) bus 
2 \ Ox; Ox, 2 \ Ox; ðr Ox3 


2 Ox; Ox, 
1 Ou, Ou) 1 ( Ous x2) 
= 一 -一 -| 一“ 一 一 一 一 0 -一 | 一 一 一 -一 一 一 
2 Ox, Ox, 2 Ox, Ox ° (2.11) 
4 (aus) —+( Qe— Su) 0 
2 Ox; Ox, 2 Ox, Ox 


如 此 分 解 的 变形 张 量 的 物理 意义 可 作 如 下 的 理解 : 反对 称 部 分 本 
来 具有 三 个 分 量 


1 (Ou; _ u) 1 ( Ou _ Ou) 1 ( Ou, _ Ou 
2 (Se ou), 2 (28 St), 2 (22 eet), 
很 容易 证 明 这 三 个 分 量 就 是 流体 微 团 绕 通 过 其 中 心 而 平行 于 三 坐 


标 轴 的 三 轴线 的 角 转 动 牵 分 量 。 显然 角 转 动 不 会 使 流体 微 困 发生 
应 变 。 所 以 所 有 的 应 变 必 由 对 称 部 分 而 生 . 已 经 证 明 


Ox1” Ox, Ox 
是 拉 应 变 ,因而 其 余 的 
4 (24 4 On), 4 (2a + us), 1 (2u + 2u) 
2 Ox Oxy 2 Ox3 Ox2 2 Ox, Ox3 


FED) EASES, SK LA ARE A EK Ou, /Ox, 那 
样 作 一 详细 计算 ， 每 一 切 应 变 具 有 两 个 对 应 的 切 应 力 ， 例 如 


2 Orz Ox, 

有 my 各 We 所 以 单位 时 间 单 位 体积 的 总 功 是 

Ou; 

ii 

Ox; 

于 是 ,最 后 得 能 量 守重 方程 
De 9 1 ôt q 1, Ou 
Dt 9 P Ox; + ep 8 加 (2.12) 


这 个 方程 可 以 用 动力 学 方程 (2.7) 改写 一 下 .。 把 每 个 动力 学 方程 
分 别 乘 以 相应 的 分 速 ;然后 加 超 来 ,得 


式 中 地 wi 当然 是 单位 质量 流体 的 动能 。 把 上 面 这 个 方程 和 方 
程 (2.12) 相 加 , 得 


P(e 十 二 wu] = Q + UiXi — 1 ða: 
Dt 2 P Ox; 
+42 (rw), (2.13) 
p Ox; 


SERA FESS CR PS BB 5S Bh BE Dt PEE, 
方程 (2.13) 的 最 后 一 项 可 以 用 方程 (2.3) 和 (2.4) 折 定义 的 粘 
性 应 力 张 量 rj 疏 写 为 
1d 
p Ox; 
Wrote ARESE E 步 演化 
D/P? 1 ð 
-petia 
所 以 能 量 方程 的 又 一 形式 是 
Op 


2 (atte eTo tenaa 
Dt Pp Ot 


1 ðq; 1 
一 一 一 和 十 一 一 (Tjuj), 2.14 
9 Ox; 3 ç 了 i) ( ) 


此 处 = e + (p/p) FAMER, WA(2.14) RAR Sa 
AATEC, WROD XE Y, Bp 


一 一 到 一 ， (2.15) 
Ox; 
那 未 能 量 方程 还 可 以 变 一 个 形式 , 这 时 方程 (2.12) 可 以 写成 


petat) o + 42e 4 8A 
Dz 2 p ôt Oz 


xiu) = L (rju) — LS (pad), 
p xi e Ox; 


1 Og: | 
一 -一 2.16 
p Əz, p Ox, -9_(r， iti). 《 ) 


如 此 形式 的 能 量 方 程 [1 ] 的 左 侧 是 烩 .动能 和 势能 之 和 的 变化 率 . 


A.3. 粘性 应 力 和 热 通 量 
在 气体 动力 学 问题 里 , 未 知 量 计 有 热力 学 参数 p, p MT, $ 


度 矢 wi, ERR g; 和 粘性 应 力 张 量 r 当 热 力学 参数 确定 时 ， 
` 内 能 。 TS 4 也 随 之 确定 .到 现在 为 止 所 有 推导 出 来 的 基本 方程 ， 
郎 连续 性 方程 .动力 学 方程 .能 量 方 程 和 流体 的 状态 方程 还 不 足以 
， 求解 问题 ， 因 为 未 知 量 的 数目 超过 现 有 方程 数目 。 和 欲 完 全 确定 问 
a, eae AE, 
由 于 从 微观 的 观点 看 来 ， 任 何 一 种 流体 都 是 由 极 多 的 分 子 所 
组 成 ， 所 以 气体 动力 学 问题 的 基本 人 处理 得 通过 统计 力学 的 原理 进 
行 ， 假定 经 典 力学 可 用 , 工 设 质点 的 疮 数 是 NN, 每 一 质点 的 自由 度 
` 数目 是 >。 这 六 个 质点 的 特定 瞬时 状态 是 为 xN 个 位 置 坐标 和 nN 
个 动量 坐标 组 成 的 吉 物 斯 (Gibbs〉 相 空间 ( 共 是 22N HE) 中 的 一 
点 所 标明 的 。 统 计 力 学 的 基本 问题 是 要 决定 NN 个 质点 的 系统 在 相 
空间 的 任何 一 点 每 一 时 到 的 几率 , 或 者 说 要 决定 相 空 间 的 2nN AK 
标 及 时 间 * LARA He, 一 旦 解决 这 一 点 ， 那 未 流体 的 任何 
一 种 宏观 性 质 ,如 热 通 量 和 应 力 , 便 都 可 以 通过 平均 几率 分 布 画 数 
RH, 最 近 欧 文 《Irving) Greve (HR (Kirkwood) 鲁 作 过 这 
个 平均 计算 . 
决定 几率 函数 的 一 般 问题 自然 是 极其 困难 的 。 好 在 处 于 常 压 
及 低压 下 的 气体 ,其 分 子 密 度 是 如 此 之 低 , 只 考虑 分 子 两 两 相遇 就 
行 了 . 换 旬 话说 , 当 两 个 分 子 接近 到 相互 有 显著 影响 的 近 距 高 时 ， 
可 以 认为 气体 中 的 其 他 分 子 离 这 两 个 相互 作用 的 分 子 十 分 造 远 ， 
因此 影响 不 到 它们 。 气 体 的 分 子 运动 论 就 是 建立 在 这 两 两 “ 相 播 ” 
的 概念 基础 上 的 .。 这 种 理论 的 几率 分 布 汞 数 还 可 做 相应 的 大 简 
化 ， 即 只 考虑 一 个 分 子 的 几 榨 分 布 范 数 ， 这 个 分 布 画 数 是 由 一 个 
名 叶 玻 耳 兹 曼 方 程 的 积分 -微分 方程 所 决定 。 恰 浦 门 (Chapman) 
ALBIS (Enskog) ([3], 第 三 章 ) 假定 分 布 画 数 随 空间 和 时 间 
的 变化 都 很 小 ,算得 热 通 量 和 粘性 应 力 张 量 是 


67 
， ， 3.1 
q k On, (3.1) 
Ou; Out) ( » 2 ) ður 
Tij = 一 一 : 十 一 一 十 -一 一 6; 3, 
、 (Se Ox, Pg ET on? (3.2) 


RHAI k ERR, 4 ERBERK, o 是 膨胀 粘性 系数 
所 以 & 和 剪 应力 及 拉 应 力 都 有 关系 ;但 往往 只 有 剪 应力 是 重要 的 ， 
因为 散 度 Ou,/Ox, 和 前 应 变相 比 往往 很 小 。 但 ye 却 只 和 拉 应 力 
BR, KERKEN, ww OBES. PMA u 是 体积 膨胀 率 或 
收缩 率 的 粘性 系数 。 A 和 jy 为 分 子 与 分 子 之 间 的 相互 作用 规律 
及 分 子 的 性 质 所 决定 的 。 例 如 , 当 分 子 没 有 内 在 自由 度 时 ,或 当 内 
在 运动 没有 激发 起 来 时 , y ER. 多 原子 的 分 子 ， 其 5 FAR. 
值得 注意 的 是 ,前 面 从 宏观 角度 论证 过 的 应 力 张 量 的 对 称 性 ,在 分 
子 运动 论 里 还 是 成 立 的 . 

这 两 个 粘性 系数 和 pe 主要 是 温度 的 函数 , 压强 对 它们 只 有 
很 微 能 的 影响 。 VAN 4 随 温度 增高 而 境 大 , 这 个 性 质 和 液体 的 
完全 不 同 。 有 关 粘 性 系数 和 热传导 系数 的 详细 讨论 可 参看 本 丛书 
的 第 IS Di. 不 过 方程 (3.1) 和 (3.2) 只 适用 于 常 压 下 的 气体 . 
在 低压 下 ， 热 通 量 和 粘性 应 力 都 出 现 新 的 项 ， 这 些 项 鲁 由 伯 尼 特 
(Burnett) ([3], 第 15 章 ) 算出 来 过 ， 十 分 复杂 。 (关于 这 个 问题 
MPs, TASS Ie, H 编 .) 

方程 (3.1) 和 (3.2) 说 明 当 温度 梯度 和 速度 梯度 为 老 时 ,不 存 
在 热 通 量 和 粘性 应 力 。 只 有 当空 间 变 化 率 存 在 时 ， 这 使 得 流体 从 
一 点 流 到 另 一 点 就 需要 活 应 新 的 情况 ， 热 通 量 和 烙 性 应 力 才 会 出 
Sh, 显然 这 样 的 活 应 并 不 可 能 瞬时 完成 ， 而 是 经 过 分 子 的 多 次 确 
捞 才 能 完成 的 。 所 以 正 是 在 流 场 中 不 同 点 企图 适应 不 同 的 热力 学 
平衡 这 样 一 个 过 程 才 产生 热 通 量 和 粘性 应 力 的 。 当 流体 状态 只 有 
很 微小 的 空间 变化 时 , 那 未 流体 就 简直 没有 调整 其 状态 的 必要 ,而 
且 很 接近 于 热力 学 平衡 ， 所 以 在 热 通 量 和 粘性 应 力 直 到 几乎 近 于 
霉 时 ， 流 体 的 状态 , 特别 是 压强 密度 和 温度 之 间 的 关系 , 必 是 热 
力学 平衡 之 下 的 关系 。 这 就 是 为 什么 应 力 张 量 r; 可 以 分 解 为 压 
aR P 和 粘性 应 力 张 量 zi 的 绿 故 ,这 个 压强 是 按 流体 处 于 热力 学 平 
衡 之 下 的 状态 方程 计算 的 . 

UHLIE (Grad) RI TPAR ALES A 
方法 .他 的 结果 适用 于 象 激 波 内 部 那样 的 条 件 ， 几 率 分 布 画 数 随 


空间 和 时 间作 迅速 变化 的 情况 .在 这 种 分 析 法 中 , 热 通 量 矢 g 和 
应 力 张 量 rr 不 能 表 为 速度 及 温度 的 导数 的 显 式 ， 另 一 方面 他 提 
出 的 附加 方程 又 包含 了 所 有 的 未 知 变量 。 所 以 在 这 种 理论 里 变数 
2 和 rz 和 其 他 的 未 知 量 p, o, T 和 u: 都 处 于 同等 地 位 . 格拉 德 的 
理论 在 稀薄 气体 动力 学 中 的 应 用 是 特别 重要 的 (第 MA, H 编 ). 
当然 分 子 运动 论 不 仅仅 给 出 方程 (3.1) 和 (3.2)。 宅 还 决定 热 
传导 系数 和 粘性 系数 与 分 子 属 性 的 关系 (第 I 千 ,D 编 )， 如 果 只 需 
要 热 通 量 和 应 力 张 量 与 其 他 指定 的 流 场 变数 之 间 的 关系 式 的 话 ， 
那 末 直接 根据 量 纲 分 析 及 坐标 变换 的 不 变性 ， 一 般 已 足以 推导 出 
来 . 例如 ,如 果 规 定 应 力 张 量 为 分 速 的 罕 间 导数 的 线性 画 数 , 那 末 
就 可 以 唯一 的 得 出 具有 未 知 参数 & 和 py' 的 方程 (3.2) KR, we AO’ 
数值 待定 ， 特 和 鲁 斯 待 尔 (Truesdell)551 鲁 作 此 元 长 计算 , 其 结果 表 
达 式 极为 复杂 .其 初始 几 项 自然 就 是 方程 (3.1) 和 (3.2) 所 给 出 的 ， 
分 子 能 量 的 三 种 主要 形式 是 平移 能 、 转 动能 和 振动 能 ， 三 者 
PIA BREAD BR CEIS. H 编 )。 所 以 如 果 气 体 状态 发 生 
疮 又 变化 , 则 振动 能 一 般 落后 于 其 平衡 什 ， 高 温 气 体 , 具 有 可 观 的 
分 子 振动 能 。 这 时 振动 能 的 滞后 是 第 二 粘性 系数 w 的 主要 来 源 ， 
不 过 这 个 问题 还 有 另 一 种 看 法 ， 这 是 堪 特 劳 维 兹 (Kantro- 
witz)11 首 先 提出 的 ， 当 激发 分 子 振动 所 需 的 时 间 远 较 激发 分 子 其 
他 形式 的 运动 所 需 时 间 为 长 时 ,这 种 方法 特别 适用 。 其 所 以 如 此 ， 
是 因为 粘性 作用 在 运动 方程 中 是 以 速度 导数 出 现 的 ， 所 以 是 局 部 
作用 ,不 可 能 表现 任何 长 时 间 的 积分 影响 ， 合 <(T) 代 表单 位 质量 
在 热力 学 温度 T 时 的 平衡 内 能 ，evin(T) 代表 单位 质量 在 T 时 振 
动能 的 平衡 态 值 ，evis 代表 实际 的 振动 能 。 假定 其 他 形式 的 内 能 
实际 处 于 平衡 状态 , 那 末 实际 的 内 能 便 是 
e=e(T)—en(T) +e, 一 eT) +Lew—ew(T)], (3.3) 
ev 一 eiT) EHAR., 现在 对 小 谐 后 量 来 说 , 振动 能 向 平衡 态 
的 趋 近 牵 可 以 近似 地 玫 为 滞后 量 的 线性 画 数 BD 
Devi 
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式 中 8 是 常数 ,具有 时 间 的 倒数 量 纲 。 把 方程 (3.3) 代 和 人 能量 方 
程 (2.12) 和 (2.13),， 从 而 就 引进 一 新 的 因 变 量 ewibp。 不 过 这 时 方程 
(3.4) 就 是 为 这 个 未 知 量 而 增添 的 一 个 方程 . 不 用 说 ， 这 样 做 之 
后 ,方程 系 中 的 4& 值 和 p 值 应 该 不 包括 振动 的 滞后 效应 在 内 了 . 
否则 同一 效应 就 算 了 两 双 .，。 所 以 w SY AR ee 
BAS, Me 的 适当 值 则 由 通常 的 分 子 运动 理论 计算 (第 I 窟 ， 
D$). testi a OSCR TE, SRA EM, Bi 
WORF Hee Rips, —- SESE Se BATE otk AS BD A AB 
达到 普通 粘性 单独 作用 时 的 两 倍 . 

在 普通 条 件 下 的 气体 ,用 方程 (3.1) 和 (3.2) 两 式 可 以 给 出 足 
够 准确 的 热 通 量 和 应 力 张 量 。 将 这 两 式 代入 动力 学 方程 和 能 量 方 
程 ,未 知 量 的 数目 恰好 等 于 方程 的 数目 ,因此 基本 方程 组 齐全 。 这 
组 方程 称 为 纳 维 埃 - 斯 道 克 斯 方程 ， 当 然 这 种 称 法 是 并 不 严谨 的 。 
本 节 以 下 讨论 都 以 纳 维 埃 - 斯 道 克 斯 方程 为 基础 。 具 体 地 说 ,这 时 
能 量 方程 (2.12) 变 成 了 


De _ 1.0(,90T\_ ? Ou; , D 
p ?Tsp Ox; (i 9) p Bx p (3.5) 
式 中 的 0 称 为 耗 散 画 数 ,其 表达 式 为 
-ai 1, (Ou Ou 
® Ta Ox; 2 ru (Ses + 2u) 
(ou buy 一 2 ay 
2 (Se t 5) + (e 3 n) (Se) (3.6) 
所 以 @ 是 粘性 耗 散 所 产生 的 热量 。 方程 (3.5) 的 另 一 种 形式 是 
Di _o+L8(06rN 1 2 
Dt oF P Ox; (A ar) + e Dt + p’ (3.7) 


式 中 4 当然 就 是 单位 质量 的 答 . 

建立 了 微分 方程 组 之 后 ， 气 体 动 力学 问题 的 完全 确立 就 只 缺 
边界 条 件 了 .边界 条 件 包 括 在 固体 物 面 上 对 速度 分 量 u ET. 
或 热 通 量 矢 g; PERRE, EDL, 流体 与 固体 物 接 触 表 面 上 
的 条 件 贫 是 个 颇 有 等 论 的 主题 ([7]，, 第 676 R). MERE TER 


体 动 力学 里 ,流体 可 以 用 分 子 运 动 论 作 微 观 处 理 . 其 结论 是 : 在 通 
常 密 度 下 ， 物 体 表 面 上 的 流体 速度 必须 与 物 面 的 速度 相等 。 he 
之 ,流体 与 物 面 之 间 的 相对 运动 必 为 零 , 此 乃 所 谓 无 滑动 条 件 ， 同 
理 , 流 体 的 温度 T 必须 等 于 物 面 的 温度 ， 这 就 是 物 面 上 无 “温度 突 
BR SRE, BRS, 当 密度 很 低 时 , 物 面 条 件 应 作 修改 (第 I, H 
编 ), 但 对 普通 密度 来 说 ， 无 滑动 和 无 温度 突 跃 是 决定 问题 的 正确 
而 又 充分 的 条 件 . 

在 上 面 的 讨论 里 ， 贫 默认 流体 可 作 单 纯 一 种 气体 来 看 ， 如 果 
流体 是 由 几 种 气体 非 均匀 地 混合 而 成 ,或 者 组 成 气体 之 闻 有 化 学 
变化 , 那 情 况 就 大 为 复杂 了 (第 工 害 ,F 编 )。 此 外 , 即使 流体 起 始 
是 均匀 混合 气 ,在 有 巨大 的 温度 梯度 的 地 区 , 热 扩 散 之 类 的 特殊 效 
应 也 会 破坏 混合 气 的 均匀 性 的 ， 好 在 飞行 器 和 流体 机 械 所 过 到 的 
主要 气体 动力 学 同 题 工 没有 如 此 复杂 的 情 驳 , 较 简单 的 纳 维 埃 -斯 
道 克 斯 方程 组 已 足够 用 . 


A.A. 基本 方程 的 积分 形式 


某 些 特 殊 问 题 用 基本 方程 的 积分 形式 比 用 前 一 节 的 微分 形式 
往往 来 得 方便 。 这 积分 区 域 是 一 固定 的 空间 域 。 设 其 体积 为 VV， 
这 个 域 可 以 是 单 连通 的 ,也 可 以 不 是 单 连 通 的 ,可 以 为 一 个 曲面 所 
包围 , 也 可 以 为 多 个 曲面 所 包围。 分 这 个 域 的 微 元 体积 为 dV, 包 
洲 的 界面 的 微 元 记 为 dA, 全 n 为 表面 微 元 dA 的 单位 法 向 矢 
(图 A.4 a)。 将 微分 形式 转换 为 积分 形式 的 基本 运算 是 引用 如 下 
所 示 的 斯 道 克 斯 定理 : 设 记 是 任 一 向 量 , 刚 


|, Sear = | Gin)a4. (4.1) 


把 方程 (4.1) 应 用 于 连续 性 方程 (2.5), 得 
0 


S | pdv + | olum) dA = 0, (4.2) 


这 个 方程 的 物理 意义 是 很 明白 的 : olun dA 是 通过 曲面 微 元 dA 
流出 这 个 域 的 质量 流量 。 所 以 对 4 的 积分 就 是 党 流出 质量 。 因 假 


TEE ea A Bl CAR, Br SA EE Fa Sah PE 
率 所 平衡 。 方程 (4.2) 就 是 表达 这 一 关系 的 。 在 一 般 有 流体 加 入 
或 抽出 的 情况 , 方程 (4.2) 的 右 侧 应 等 于 此 域 中 的 源 和 汇 的 强度 的 
代数 和 。 


图 A.4 a 


同 理 ,动力 学 方程 (2.6) 转 换 成 为 


2| pu; dV + | pu;(ujn;) dA 
Or Jv A 


= | pX;dV + \ (Cajn;) dA, (4.3) 
f, | 


这 个 方程 的 物理 意义 也 是 很 简单 的 。 方程 的 左 侧 是 此 域内 的 流体 
在 x; 方向 的 动量 增 量 ， 而 右 侧 则 是 产生 这 个 增 量 的 原因 , 6k 
力 的 合力 和 表面 应 力 之 合力 。 方程 (4.3) 当然 是 基于 VV 域内 没有 
加 入 流体 ,也 没有 抽出 流体 的 这 个 假设 的 。 当 此 假设 不 成 立时 , 方 
程 (4.3) 的 右 侧 还 得 加 上 由 于 引进 或 抽出 流体 所 引起 的 动量 变化 
Æ, 

除了 线 动 量 之 外 ,还 可 以 计算 对 三 条 轴线 的 角 动 量 。 为 此 ,最 
好 引用 一 个 符号 6ix。， 这 个 符号 的 意义 是 : 

Sin = l, AIÈ i AGAR, Mi, 7, RK 1, 2, 3 的 顺序 轮转 ; 

Sak = —1, WR i AGAR, Mi, 7, k IER 1, 2, 3 的 顺序 输 

转 ; 


© 16 。 


wi bet 


Sik = 0, 其 余 情 况 . 
例如 ôi = 1, 而 51% = 一 1， 方 程 (2.6) 可 以 号 成 


pur 0 — Onn: 
Cour 4 2 D) = pX, + CH, 
Ərz Oz, (purui) OAR Ox; 


以 Opjk Xi 乘 此 方程 ， SRB ARR VERS, 


2j põik Xiur dV + | pOik xjurlum) dA 
Dr Jv A 


= |, oda xjX,aV 十 |, Oak x(t) dA, (4.4) 

这 是 绕 x; 轴 的 角 动 量 方 程 ， 式 中 各 项 依 取 和 约定 对 i MR 取 和 . 
RAGS th, ARIA V EAV pih, 方程 (4.4) 还 必须 泛 
当地 包括 这 些 源 和 汇 的 影响 , 

为 了 把 能 量 方 程 (2.13) 变换 成 积分 形式 , 可 以 看 到 利用 连续 
方程 (2.5) 可 将 方程 (2.13) 写 成 : 

1 0 1 

-— p (e + E umi) + Dn, OM (e + > uju) 


Xi 
04: 
= 00 + puiXi — SE + 9 (oy), 
IAM? V ASAT AE FF 
| 
By vy? e + 2 Uit; az 十 | p “十 了 ujt; | un; dA „PEN 


+ | ou;X;dV 一 | (gin;)dA + | (ziuini) dA, (4.5) 
v A A 


Ze MURR EAR AY AN BE FSB Aes, 右 侧 表示 这 一 增 量 
的 来 由 。 第 一 项 是 加 入 的 热量 ,如 通过 化 学 反应 ,第 二 项 是 彻 体力 
所 做 的 功 , 第 三 项 是 通过 传导 所 失去 的 热量 ,第 四 项 是 表面 4 上 的 
应 力 所 做 的 功 。 能 量 方程 男 一 种 积分 形式 可 以 得 自 方程 (2.16)， 
gp 


2f ( toe ) 
Br vy? h+ its +H dV 
+ | p(s ttum ) uma = | oQ aV 
4 2 v 
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+2 pav +2] oy dV -| (qmi) dA 
Or Jv Oz Jv A 


+ | (Tauni) dA, (4.6) 


FBAWMRARIG. BPESA NLRA, 右 侧 表 示 这 个 增 
量 的 来 由 ， 第 一 项 仍 是 加 入 的 热量 9, 第 二 项 是 压力 能 的 变化 这 ， 
第 三 项 是 势能 的 变化 牵 , 第 四 项 是 通过 传导 所 损失 的 热量 损失 率 ， 
最 后 ,第 五 项 是 表面 4 上 粘性 应 力 所 做 的 功率 . 

对 积分 形式 的 基本 方程 可 以 作 如 上 的 简单 理解 ,事实 上 ， 如 
果 愿 意 的 话 ,可 利用 积分 形式 的 方程 作为 讨论 的 起 点 ,再 从 该 方程 
导出 微分 方程 . 

绕 物体 的 流动 . 以 流体 绕 一 物体 的 定常 流动 (图 图 A.4b) 为 例 
来 说 明 积 分 方程 的 应 用 。 ERAAN 4 与 物体 闻 的 空间 为 流体 
域 , 曲面 4 包围 物体 在 内 这 样 , 流体 域 的 表面 4 便 有 两 部 分 : 


图 A.4b 

Bl 4, 和 物体 表面 4;， 进 一 步 假 定 没有 彻 体力 , 且 9 = 0; 这 时 动 
量 方程 (4.3) AAS TARE., 由 于 物体 表面 上 流体 速度 
AS, 所 以 方程 (4.3) 的 左 侧 第 二 个 积分 中 沿 A 的 动量 的 通 量 积 
PAE, woo, 根据 作用 力 与 反作用 力 相 等 的 原理 , 知 方程 (4.3) 
的 右 侧 第 二 个 沿 4. 的 积分 应 等 于 作用 在 物体 上 的 合力 Fo 但 与 
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它 方 向 相反 ， 所 以 
| ou;(ujn;) dA = | (min) dA Fy, (4.7) 


IRER, BTA KP E 41 的 积分 求 得 物体 所 受 的 
作用 力 F. 一 个 好 用 的 办 法 是 取 曲 面 41 远离 该 物体 ， 以 使 上 列 
积分 计算 大 大 简化 . 

用 同样 的 办 法 于 角 动 量 各 方程 (方程 4.4), 可 得 作用 于 固体 物 
的 绕 坐 标 轴 的 力矩 M; 的 方程 为 


| OÕijktjul umn) dA = |, Skr; rrini) dA — Mj, (4.8) 
记 固 体 物 所 吸收 的 热量 为 互 , 利用 方程 (4.6) 则 得 能 量 方 程 如 下 : 
L p (z + 7 wu) ujn, dA 
=—H+| Cwmda—| Camda, (49) 
如 用 粘性 应 力 rz 则 方程 (4.7),(4.8) 和 (4.9) 可 以 写成 


F; = |, [Tun 一 pni 一 pui( uin; )] dA, (4.10) 
1 
i7 |, Bigxi [Tini 一 png 一 ou,Cuyn,) | dA, (4.11) 


H = | [eam — Qin; 一 0 (z + + wu) “| dA. (4.12) 
Al 


这 三 个 方程 表明 :作用 在 物体 上 的 力 就 是 沿 A1 的 应 力 积分 和 治 
A 的 动量 积分 之 差 ; 同样 的 ,作用 在 物体 上 的 力 挎 等 于 沿 41 应 力 
的 力矩 积分 和 动量 矩 积分 之 差 ; 最 后 , 物体 所 吸收 的 热量 等 于 4 
面 上 的 粘性 应 力 所 做 的 功 与 热量 、 烩 及 动能 三 者 流出 界面 4, 的 总 
外 流量 之 差 . 


A.5. 相似 性 和 流动 参数 
为 清晰 起 见 , 区 考虑 一 具体 问题 . 设 在 未 扰 速 度 为 UU 并 具有 


1) 另 一 种 讨论 相似 性 的 方法 见 第 VIL, D. 


重力 作用 的 无 限 流 场 里 , 存在 一 典型 的 线性 尺寸 为 工 的 物体 ,， 工 
可 以 是 物体 的 长 度 . 物体 正在 作 频 率 为 了 的 振动 。 自由 流 中 的 各 
物理 量 注 以 下 标 co , 无 量 纲 的 量 注 以 星 号 ,于 是 

l ox 


一 了 T = TaT*, X: = eX}, 
x; = Lef, wi = Uwur, tj 一 tell h, 
7 (5.1) 
0 = pop", e= eue", di 一 个， 
L 
U oo 
P= peP*, O= —=O*, 
此 处 
* k OT* 
i= 5.2 
q he Ox* > ¢ ) 
* * , * 
th —+ (2t + our) + (4-2) sy Ott, (5.3) 
Hæ Ox} Ox} Loo 3 Ho Ox*¥ 


8 是 重力 常数 .利用 这 些 无 量 纲 变数 , 则 连续 方程 (2.5) 可 以 写 为 


Lf\Op* | Op*ur _ 
(Gt) Se + ee o, (5:4) 
动力 学 方程 (2.9) 变 成 
(H) 4. yp Set. (_to.) 4 20" 
D/O Oxf Poll?/ po* Ox* 
+ (25) X7 + (==) A Ori (5.5) 
U? PoU L/ p* Ox¥* 


记 定 压 比 热 为 cp, 其 在 自由 流 中 的 值 为 (c*)-。 刚 能 量 方程 (2.12) 
成 为 


v/o at 
一 Or 一 (一 ) (<=) (eTa) 1 Oa 
Cpeoftc4 \PoolJ L ĉo / p* Ox 


— (Pe U? 所 Out + (12 u? tH Out 
2 * * . (5.6) 
PoU? Sea p* Ox} PoU L/ \ew/) po* Ox* 
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LAA AEE A RFR, CER EP IE 
问题 的 参数 。 其 中 L/U BPR “RARE”, U2/ gL AMHR 
数 ,而 o-U 工 /we 为 雷诺 数 ， 比 值 各 /pu 具有 速度 平方 的 量 纲 ， 事 
实 上 NV pool Peo 是 和 分 子 无 规则 运动 的 平均 速度 同一 数量 级 ， 也 和 
自由 流 中 扰动 传播 速度 同一 量 级 。 所 以 U/W pe/pe 是 和 自由 流速 
与 扰动 传播 速度 (TR) 之 比 同一 量 级 。 当 然 流速 对 声速 之 比 ， 
就 是 马赫 数 。 所 以 po/pwU0? 就 是 马赫 数 平方 的 倒数 。 eo 和 分 子 
骚动 的 动能 同 量 级 . 因而 eo 双 和 分 子 的 平均 速度 之 平方 或 声速 
的 平方 同 量 级 . 于 是 U2/ew 又 是 马赫 数 的 平方 。 因子 poe po/he 
是 传 热 效应 和 粘性 效应 的 相对 重要 性 之 度量 ， 这 个 数 称 为 普 翩 伦 
数 . 因数 cywTw/ew WETS ARELI, EBON 1 的 一 个 数 
值 . 
然而 流体 力学 的 一 个 非常 重要 的 问题 就 是 相似 问题 ， 具 体 的 
说 , 我 们 所 关心 的 是 一 些 几 何 相 似 的 物体 , 相似 物体 可 以 用 平移 ， 
或 毛 比 例 放 大 或 缩小 的 办 法 使 之 彼此 重合 。 再 考虑 前 面 所 壕 在 无 
限 流 场 中 的 物体 振动 问题 。 问题 是 :对 于 两 个 几何 相似 的 物体 应 
该 如 何 规 定 其 自由 流 参 数 UU，p。， pos Tw, Ho, ko, co A L 5f, 
才能 使 两 组 参数 都 导出 同样 的 无 量 纲 微分 方程 ( 即 方 程 5.4, 5.5 
和 5.6) 来 9 如 果 能 解决 这 个 问题 的 话 , 那 未 对 应 于 这 两 组 参数 的 
两 个 流 场 就 是 “相似 ”的 ， 也 就 是 说 ,对 于 两 流 场 说 来 ,如 把 一 切 变 
数 都 改 成 无 量 纲 的 话 ， 变 数 之 间 的 画 数 关系 两 者 都 是 一 样 的 ， 这 
样 就 使 气体 动力 学 问题 得 到 大 大 的 简化 ， 可 以 节省 许多 时 间 和 精 
力 。 这 样 的 收获 自然 不 可 能 不 付出 一 点 代价 ， 那 就 是 得 辆 牲 一 些 
普通 性 ， 换 句 话说 ,不 对 流体 的 属性 加 上 某 些 限制 的 话 ,这样 有 力 
的 相似 律 是 得 不 出 来 的 。 好 在 在 常温 和 常 压 下 气体 足 可 以 近似 地 
当 作 “ 完 全 气体 "处理. 完全 气体 是 遵循 下 烈 状态 方程 的 气体 : 
P = RoT, (5.7) 
式 中 的 R 是 特定 气体 的 气体 常数 ， 等 容 比 热 c, 和 等 压 比 热 cs 都 
可 能 是 温度 的 画 数 ,但 二 者 之 间 存 在 下 型 的 关系 
cp =ot+ &, (5.8) 
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MR y 是 两 种 比 热 之 比 ， 


S 


y=, 
Cy 
HRED.) E EREA 
pay 81. 
于 是 内 能 。 成 为 
T 
e = f c dT, 
于 是 
和 
C pol œ 0 € p00 Yoo 
及 
ee 
cm Cee) 
如 果 M 代表 马赫 数 , 那 末 
eel = yuM2 
及 


2 
ut = (Yæ — 1)M2, C poo Teo 


oo Co 


(5.9) 


(5.10) 


(5.11) 


(5.12) 


(5.13) 


记 雷 诺 数 为 Re, BERN Pr, ERER H Fr, tRYRMAEA 上， 


ep 
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(5.14) 
(5.15) 
(5.16) 


(5.17) 


则 连续 方程 ,动力 学 方程 和 能 量 方程 分 别 变 成 
Bo” + Op" u¥ 


=0 5.18) 
“ort Ox? ° € 
p Ove +u x uf _ 1 16 1 Lys 
* Oye “i Ox} YoM?, p* Ox Fr “ 
1 OT; ( 
i 5.19) 
Resp” Ox¥” 
k Oe* 十 ut ðe* _ yw 1 (se ) E Ogr* 
Oz* , Ox* PreoRex Co e* Ox} 
_ Ye al (ete) t Our 
Yoo Coo p* Ox} 
十 (yw 一 LMS ( eer) TÀ Do (5.20) 
Rew ĉo / p* Ox¥ 
状态 方程 (5.7) 可 以 写成 如 下 的 无 量 纲 形式 
p* = o*T*, (5.21) 


研究 一 下 (5.18) 到 (5.20) 各 方程 以 及 (5.12) 和 (5.13) 各 方程 ， 
可 以 看 出 来 , 欲 使 两 流动 相似 , 其 参数 值 必须 满足 下 列 要 求 , Blam 
GEE e, 马赫 数 Mo, BiG Reo 和 傅 劳 德 数 Fr 相同 。 此 外 ， 
流体 的 属性 还 必须 使 两 种 流动 的 普 良 伦 数 Pr 和 两 种 比 热 之 比 7 
都 相同 。 还 有 , 比 热 , 热传导 系数 & 和 粘性 系数 上 和 py' 随 温度 的 
变化 规律 都 必须 使 cp/cpw, R/ Roos tf poo 及 p/pw 各 比值 各 成 为 温 
BET* ME, 如 果 cp, ko uo u 随 温度 了 的 霜 变 化 的 话 , 最 
后 这 个 条 件 便 能 得 到 满足 。 不 用 多 说 , 方程 (5.1) 是 包含 了 这 样 一 
条 假设 的 , 部 加 给 单位 质量 的 热能 2 是 和 Te。/ 工 或 IT。/ 工 成 正 
比 的 。 当 上 列 的 这 些 条 件 全 部 满足 的 时 候 ， 那 未 这 两 种 流动 就 有 
相同 的 无 量 纲 微分 方程 组 .。 只 要 二 者 的 无 量 纲 边界 条 件 也 相同 ， 
那 末 二 者 用 无 量 纲 变 数 表 达 的 解 应 是 一 样 的 . 

重要 的 边界 条 件 是 物体 表面 的 条 件 。 这 就 是 流体 和 固体 物 面 
之 间 相 对 速度 必 为 零 的 条 件 ， 以 及 物 面 上 的 温度 条 件 。 流速 的 条 
件 是 一 种 齐 次 的 条 件 ， 它 不 会 再 引出 新 的 相似 参数 来 ， 但 物体 彭 
面 的 温度 是 规定 物体 与 流体 之 间 的 热 通 量 的 ， 将 引出 一 个 新 的 相 
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似 条 件 来 ， 合 (gi)w 代表 物 面 上 的 热 通 景 矢 .。 则 这 个 边界 条 件 的 
无 量 纲 形式 应 是 


Qu = ee -gw pry) ys t"), (5.22) 
式 中 的 oy 是 物 面 各 点 的 无 量 纲 坐 标 、 要 相似 , Eee FS AE 


fl, 根据 这 一 个 条 件 ， 壁 面 无 量 纲 温 度 应 确定 为 T 名 一 Ty/T。. 
于 是 条 件 (方程 5.22 ) 就 等 于 


一 ¢ Pw 一 Cgi)wL = x? * 
Nu T31 RT, ~TO Nu(Ca¥)w, t*), (5.23) 
Nu 称 为 努 赛 尔 传 热 数 . 


总 括 起 来 ， 对 于 一 种 非 定 常 的 、 有 热传导 的 和 粘性 的 完全 气 
体 的 运动 来 说 ， 若 单位 时 间 加 给 单位 质量 的 热量 为 0， 且 受 重 力 
影响 ， 则 相似 条 件 就 是 要 下 烈 各 参数 相同 ; 马赫 数 M。， 雷 诺 数 
Reo, toa x, SFB Fr, MBER Pr。， 两 种 比 热 之 比 
Yos 作为 T* 的 画 数 的 比值 cp/cpos k/ko, Ut/ Hos H / tos 作为 x7 F 
六 BOER BA OL/UT o 以 及 最 后 作为 (x)w 和 e 的 画 数 的 努 赛 尔 
Ze. 不 用 说 , 当 物 体 运 动 不 是 纯 振 动 的 时 候 , 只 要 以 1/r 代替 频率 
fp 上 面 的 分 析 仍 可 适用 , T 是 问题 里 的 特征 时 间 . Sc ees aS 
变 成 了 L/Ur, 如 果 粘 性 和 热传导 的 作用 都 可 以 略 去 不 计 的 话 ， 
那 末 上 述 的 相似 条 件 中 可 以 不 要 雷诺 数 、 普 朗 从 数 ,， /ko, Lk/ ie， 
/bw 以 及 努 赛 尔 数 .。 如 果 没 有 外 在 的 力 场 作用 ， 那 末 侍 劳 德 数 
RSS. URGE HN, MABE, 如果 没 有 
热量 加 入 , OMMETRB, 如 果 比 热 是 常数 ， 那 末 比 值 cp/cpw 又 
可 以 不 列 入 相 亿 条件。 所 以 如 果 要 研究 的 只 是 有 粘性 导热 的 完全 
气体 的 定常 运动 , 且 气 体 的 比 热 是 常数 ,热传导 系数 和 粘性 系数 都 
是 常数 ,工交 有 热量 加 入 , 也 没有 外 力 场 的 作用 , 那 末 流 动 的 相似 
性 便 为 马赫 数 、 雷 诺 数 、 普 朗 伦 数 、 努 赛 尔 数 以 及 两 种 比 热 之 比 所 


规定 . 
A.6. 理想 气体 


对 于 任何 的 实际 气体 或 数 种 气体 的 混合 物 ，& 和 上 的 相对 值 
。24 。 


AEEA E 1 的 量 级 ， 事 实 上 , 按 分 子 运 动 论 , 单 原子 气体 
AOR RITE Bee SET 1 的 。 所 以 粘性 的 影响 和 热传导 的 影响 同等 
BR: 作 分 析 时 考虑 了 一 种 影响 ,同时 也 得 考虑 另 一 种 影响 ,这 样 
才 是 正确 的 。 粘性 系数 上 和 u 以 及 热传导 系数 的 数值 都 很 小 . 
当 物 体 的 尺寸 为 器 的 量 级 ,速度 口 为 每 秒 数 十 中 或 数 百 中 的 量 级 ， 
压强 po 为 1 个 大 气压 左右 时 , 雷诺 数 是 十 分 巨大 的 , 约 在 一 百 万 
左右 .方程 (5.19) 和 (5.20) 说 明 一 般 情况 可 以 忽略 粘性 和 热传导 
作用 .只 在 温度 梯度 和 速度 梯度 很 大 的 地 区 ， 粘 性 和 热传导 才 会 
起 重要 作用 . 什么 地 方 的 温度 梯度 和 速度 梯度 很 大 呢 ? 在 亚 声速 
流动 里 只 有 边界 层 里 才 有 这 样 的 大 梯度 ;在 超声 速 流动 里 , 激 波 里 
面 也 有 很 大 的 温度 梯度 和 速度 梯度 。 所 以 只 有 在 边界 层 或 激 波 里 
才 需 要 考虑 粘性 和 热传导 的 作用 。 可 以 认为 在 这 些 区 域 之 外 的 流 
体 是 无 粘性 的 和 不 传 热 的 。 这样 的 流体 称 为 理想 可 压缩 气体 ， 如 
果 流 体 又 是 完全 气体 , 那 末 就 称 为 理想 完全 气体 ， 用 于 理想 气体 ， 
运动 方程 由 于 没有 粘性 应 力 和 热 通 量 项 ,而 得 以 大 为 简化 . 

在 边界 层 和 激 波 里 ， 粘 性 和 热传导 项 都 是 很 重要 的 ， 所 以 处 
理 这 类 问题 是 相当 复杂 的 (参看 第 贰 稚 D BRBIVE). 如 果 
边界 层 和 激 波 区 中 的 粘性 方程 的 解 必得 和 这 些 地 区 以 外 的 非 粘 性 
方程 的 解 连 在 一 起 ， 才 能 使 整个 问题 同时 得 到 解决 的 话 ， 那 未 序 
使 雷诺 数 很 大 ， 气 体 动力 学 问题 也 不 会 得 到 很 大 程度 的 简化 。 幸 
好 大 多 数 问 题 都 不 是 这 样 的 情况 .通常 激 波 是 如 此 之 薄 ， 可 以 认 
为 它 的 厚度 是 无 限 小 量 , 而 在 激 波 两 侧 的 气体 , 由 于 在 激 波 之 外 ， 
所 以 都 是 理想 气体 这 时 激 波 不 过 是 一 个 数学 上 的 不 连续 面 而 
B. 一旦 给 定 激 波 前 的 条 件 和 波 的 形状 , 突 跃 量 是 很 容易 计算 的 ， 
无 需 用 到 粘性 流 微 分 方程 。 这 也 就 是 说 ， 可 以 容易 地 得 到 解 的 一 
个 重要 部 分 。 边 界 层 的 问题 则 困难 得 多 ,因为 它 的 厚度 虽 很 薄 , 却 
又 不 可 忽略. Pit, TERA EMRE RE, ORT RA 
Prva SS ot AVE Bye A A A, TAREE ACER Bh AALS — Ah 
合 平 逻 辑 的 方法 是 , 先 略 去 粘性 和 热传导 ,但 如 果 需 要 的 话 必 须 包 
括 激 波 的 突 跃 ， 把 流 场 计算 由 来 ， 然 后 用 这 样 初次 近似 求 得 的 物 
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AS LHER, HEME HAKAN DFR A & 
出 边界 层 之 后 ， 就 可 以 算出 物 面 上 的 切 应 力 和 热 通 量 ， 这 样 的 迭 
代 程 序 使 问题 大 大 简化 ,这 就 是 经 典 的 普遍 伦 边 界 层 理论 的 基础 

对 亚 声速 流 来 说 上 面 所 述 分 析 方 法 是 很 成 功 的 。 在 超声 速 流 
中 ,越过 激 波 的 强大 压强 升 高 是 可 能 大 大 影响 边界 层 的 发 展 的 ,而 
改变 后 的 边界 层 双 可 能 显著 地 影响 激 波 的 形状 。 这 种 相互 影响 称 
为 边界 层 和 激 波 的 相互 作用 ， 这 是 近来 大 力 研 究 的 一 个 课题 ， 在 
高 超声 速 流 , 郎 马赫 数 十 分 大 的 流动 中 , KBR PHL, 
激 波 和 边界 层 很 靠近 ， 即 使 边界 层 很 薄 也 会 影响 到 激 波 的 形状 ， 
这 种 情况 可 以 称 为 高 起 声速 流 的 边界 层 - 激 波 相互 作用 . 有 这 类 
的 相互 作用 出 现时 ,经 典 的 边界 层 理 论 就 不 再 适用 了 . 

另外 一 种 简化 通常 总 是 可 能 的 。 倩 劳 德 数 Fr 是 /gL, 此 外 
工 是 物体 的 典型 尺寸 ,UU 是 典型 速度 ,如 果 外 力 场 是 重力 的 话 , £ 
是 重力 常数 。 所 以 飞行 器 的 气体 动力 学 问题 中 ， 依 劳 德 数 总 是 十 
分 巨大 的 ,大 到 使 动力 学 方程 (方程 5.19) 的 外 力 项 小 到 可 以 略 去 . 
这 时 外 力 就 可 以 略 去 不 计 . 

从 以 上 几 节 的 讨论 可 以 看 清楚 ， 在 气体 动力 学 中 重要 的 问题 
就 是 没有 和 外力 场 作用 的 理想 气体 的 流动 问题 ， 其 动力 学 方程 没有 
粘性 应 力 和 和 外力, 其 能 量 方程 没有 热 通 量 项 和 粘性 耗 散 项 ,因而 问 
题 就 大 大 简化 了 .和 微分 方程 的 简化 同时 ， 固 体 物 面 上 的 边界 条 
件 也 得 到 放松 。 现 在 只 要 求 相对 物体 表面 的 流速 和 物 面 相 切 ， 而 
不 需要 为 零 了 ， 物 面 上 的 温度 条 件 也 不 要 了 。. 现在 流体 的 温度 为 
其 他 的 变量 毛 单 一 决定 ,七 在 物 面 上 的 值 是 不 能 任意 规定 的 . 


A.7- 理想 气体 的 非 绝热 流动 ， 环 量 与 涡 量 


没有 外 力作 用 的 理想 气体 的 方程 组 是 


ao + OC pui) = 


Oz Ox; ? (7.1) 


-一 一 了 六 (7.2) 


D (2 十 工 ua) = 0 十 1 op (7.3) 
Dt 2 p Ot 


这 个 方程 组 货 是 希 克 斯 Hicks) PRANAJA Th, AE 
绝热 流 是 和 无 热 加 入 ( 即 9 一 0) 的 绝热 流 相反 的 一 种 流动 。 希 克 
斯 要 寻求 的 是 燃烧 所 造成 的 加 热效应 ， 所 以 他 用 非 绝 热流 作为 燃 
烧 场 中 进行 反应 中 的 气体 流动 的 模型 ， 当 运动 是 定常 的 时 候 ， 朗 
当 流 场 参 数 不 是 t MEHR, Ml 8p/6: 一 0, 方程 (7.3) 表明 对 
FRAT it Rae A WR ker ABBE AN SE CS SET TIA O. 
在 这 个 意义 上 , EREZA ATE RAS “RABEL” 2°, 
是 个 十 分 有 用 的 概念 . 在 特殊 的 绝热 流动 下 ，9 一 0， 总 能 量 沿 
任何 一 微 田 的 迹 线 ,或 说 沿 任何 一 条 流 线 , 必 是 常数 .不 过 应 该 注 
意 上 述 有 关 这 条 定理 的 限制 条 件 ， 这 结论 只 在 无 粘 不 传 热 的 气体 
作 定 常 运动 时 才 成 立 . 

男 一 个 侧 未 讨论 过 的 重要 的 流 场 变 数 是 涡 矢 岛 。 涡 矢 表 示 流 


体 微 团 的 内 旋转 的 给 度 ,定义 是 
Q=Vvxu, (7.4) 
这 里 采用 习 用 向 量 符号 , u ERR, EE u, 
u = iu, 十 ju + ku, (7.5) 


V X 是 旋 度 算 子 , 故 


Q=V7 x ui (Se — en) + j( Se — 2) 
Ox2 Oxs Ox; Ox, 


Ou, _ Ou 
+ k ( Ox, a), (7-6) 
这 些 方 程 里 的 i, j,k 分别 是 x1, xz, xs 三 个 方向 的 单位 向 量 。 拿 
方程 (7.6) 中 诸 量 和 方程 (2.11) 中 诸 量 对 比 ,可 以 看 到 , 涡 矢 的 三 个 
分 量 是 等 于 变形 张 量 的 反对 称 部 分 的 三 个 量 的 二 倍 的 ， 所 以 涡 舌 
等 于 流体 微 团 的 角速度 矢 的 二 倍 . 
用 向 量 符号 ,连续 性 方程 成 为 
Se +7 (pon) = 0, (7.7) 
ATV: 是 散 度 算 子 , 故 


. 27 。 


Y: (pu) = Sows) (7.8) 


于 是 动力 学 方程 变 成 
Du _ 4 (ye vu =—+ vp, (7.9) 
Dt Or p 
式 中 的 Vp 是 压强 ”的 梯度 矢 。 能量 方程 可 以 写成 
Di 0 44 2e, (7.10) 
Dt p Of 
式 中 的 如 是 总 能 量 ,或 
加 一 大 十 > uiii, (7.11) 


现在 根据 向 量 分 析 可 以 把 这 些 基 本 方程 作 一 些 变换 : 


v (Tum) = (Tu) = (a. vu tux (v x u). 


2 
(7.12) 
所 以 方程 (7.9) 可 以 写成 
Du 1 _ —__ ti 
oe 十 了 了 (+ ui) ux Q vp. (7.13) 
如 果 * RRA RR SB EER 
T ds = dh — = dp (7.14) 
及 
Ds — 
T Dr . (7.15) 
比 热 为 常数 的 理想 气 休 的 状态 方程 可 以 用 s, p，p RRK 
p = const pe", (7.16) 
用 方程 (7.14) 消 去 压强 梯度 项 ,动力 学 方程 就 还 可 以 写 为 
Du Tv, 
> TVs — Vh (7.17) 
或 
M ux Q = TVs — VW, (7.18) 
t 


此 关系 式 通常 称 为 克 罗 柯 (Croco) 定理 (参看 [9]).。 用 方程 
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(7.15) RE 0, 能量 方程 成 为 
Dt Dt p ôt 
利用 这 些 变 换 方 程 ,可 以 研究 理想 气体 流动 中 涡 量 的 产生 了 . 
就 方程 (7.13) 两 侧 取 旋 度 ,并 注意 VV .YY x u 一 0, 得 


8B + (a VR — (R - vu + QC - u) 


(7.19) 


=-y x (+ve) 
p 


或 
DQ _ . _ . 1 
F (Q - vu Q(V-wW-V xX (4 ve); (7.20) 
这 个 方程 的 又 一 种 形式 可 用 方程 (7.14) 推 得 
2a = (Q -Yj —Q(Y -u) + YT X Ys. (721) 


如 果 p RE 2 的 画 数 , 那 末 方程 (7.20) 右 侧 最 后 一 项 可 写成 一 个 
标量 画 数 的 梯度 之 旋 度 ， 因 而 该 项 为 零 ， 如 果 单 位 质量 的 炳 值 是 
常数 , BARAT Vs 一 0, 而 方程 (7.21) 的 最 后 一 项 为 零 ， 只 要 气 
体 是 均 炳 的 ,部 气体 到 处 的 炉 都 根 同 ,这 两 个 条 件 都 彼 满 足 ， 因 此 
在 理想 气体 流 中 , 流体 微 田 的 涡 量 的 变化 紊 是 由 方程 (7.20) 或 
(7.21) 右 侧 的 头 两 项 组 成 的 。 如 果 流 场 上 每 一 点 的 熏 在 某 一 瞬 
AA, BK D/D: 一 0， 所 以 流 场 上 每 一 点 的 涡 量 保持 为 办， 
这 种 流动 称 为 无 旋 流 ， 另 一 方面 如 果 流 动 不 是 均 糯 的 ， 朗 流 场 上 
SAS ARABI AOS, BRED EEE — BRIE], LA Q 
RER, 方程 (7.21) 右 侧 最 后 一 项 也 会 使 下 一 瞬间 涡 量 @ 不 为 需 
的 。 所 以 非 均 精 流 不 可 能 是 无 旋 流 . 因此 , 无 旋 意 味 着 均 燃 ， 然 
而 , 均 焕 当然 并 不 意味 着 无 旋 . 

现在 来 考虑 绝热 的 连续 流动 , 设 该 流动 是 由 均匀 更 下 状况 (在 
t 一 0 时 , u = 0) 开始 流动 起 来 的 。 按 方程 (7.15 ) 看 ,只 要 流动 是 
绝热 和 连续 的 ， 任 何 一 个 流体 微 田 的 业 值 都 是 常数 。 但 在 上 一 0 
时 ,到 处 炳 值 都 是 相 辐 的 ,也 邵 初 始 状 观 是 均一 的 。 因 此 运动 将 永 
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EIA, IH, 在 1 二 0 时， 因 各 处 的 流速 是 零 , He Q= 0, 
所 以 在 理想 气体 作 绝 热 连 续 流 动 的 假设 下 ， 从 议 赴 状况 起 始 的 任 
何 一 种 运动 都 是 无 旋 的 。 不 是 因为 在 超声 速 流 中 要 出 现 激 波 的 
话 , 那 这 里 面 简直 就 包括 了 所 有 的 势 流 问 题 。 在 理想 气体 里 , 激 波 
是 突 跃 面 , 它 可 以 产生 涡 ( 贿 看 第 亚 餐 , DAMER) EMEA 
激 波 存在 , 寂 所 产生 的 涡 量 通 常 是 不 大 的 ;所 以 尽管 无 旋 流 是 很 特 
殊 的 流动 ,但 却 具 有 很 大 的 实用 重要 性 . 
方程 (7.20) 或 (7.21) 右 侧 第 一 项 (Q. V)u 可 以 看 作 是 沿 涡 
矢 信 方向 的 线 或 涡 线 的 弯曲 效应 ， 为 说 明 这 一 点 , 假定 在 上 瞬间 
涡 舌 只 在 x 轴 的 方向 有 一 个 分 量 1, MAA a, 那 末 (Q -v)u 
有 两 个 分 量 
Ou Ou; 
0 FQ, Oxy 
一 个 是 方向 , 一 个 是 方向 。 现在 来 看 一 微 元 线段 (图 A.7)， 
这 微 元 段 在 t 时 和 0 重合 , EKEK da. 经 dt 时 间 后 , 线段 的 
Feta x1, x25 zx3 移 到 了 xj, x 十 wz dt, xa 十 Ws dit。 AM xit dxs 
X25 X3 移 到 了 
Ou 


xı + ax, xı 十 (un + 一 一 za dt, x3 + (us + Qe gx, ) dt, 
Oxy Ox, 


所 以 原来 的 线段 现在 变 得 和 xi 轴 相 倾 笠 了 , 它 现 在 和 x 轴 的 方向 
RL EF (Ouz,/Ox;) dt, 和 x3 轴 的 方向 余弦 等 于 (86us/8x1) dt, H 
在 假设 涡 量 0O 原来 和 x 轴 重 合 , 并 和 流体 一 起 运动 ,而 且 保 持 写 
的 强度 不 变 、 那 未 经 dt 时 间 后 , ERAT x 轴 方 向 的 分 量 , 其 值 
等 于 Q(Bw/6x1) dt, 还 有 xs 轴 方向 的 分 量 (Ous/6x1) dt, ARYA 
沿 a 轴 涡 随 时 间 的 增加 率 是 O,(0u,/On), W xs AES ae 
IMA O,(Ous/Ox), 3XTE(Q+> 了)u， 所 以 这 个 量 表 示 由 于 涡 
线 旁 曲 而 生 的 涡 增加 率 . 同 理 ,第 二 项 , 朗 和 散 度 V - u 或 体积 变 
化 率 有 关 的 一 (7 . a) 可 以 理解 为 包围 涡 线 的 流 管 的 伸延 效应 
所 以 如 果 压 强 只 是 密度 的 画 数 ,或 者 如 果 炳 是 个 常数 ,因而 流体 微 
Mis Q 的 变化 只 有 上 述 两 项 的 话 , 那 未 上述 的 结果 可 以 认为 是 渴 
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HORTA RS ia ak QO aS HOFER RSE AE, 对 空间 和 时 间 说 来 
都 不 变 。 FARA Sat, oR EUR EE BE RE. 
这 种 不 变性 说 明 涡 的 概念 是 很 有 用 的 ， 


XI, X2, X3 xı 十 dx1, x2, X3 


xt 十 dxi, 


x2 + (= + Ba, Net, 
Ox1 


iy x+ (» + a 


Ee dxdt 


x3 


x, x2 十 Uadt, xa 十 usdt 


图 A.7 


海 姆 填 尔 兹 定理 的 内 容 还 可 以 用 另 一 种 办 法 证 明 ， 这 种 证 法 
也 许 更 直接 些 : 合 T 为 沿 任 何 闭 合 路 绕 C 的 环 量 ,定义 为 


r= $ u - dl, (7.22) 


式 中 的 dl 是 路 线 C 上 的 微 元 线段 。 引用 围 线 积分 的 斯 道 克 斯 定 
理 则 


=| (xu) A=] Q- aa, (7.23) 
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式 中 的 4 是 单 连 曲面 , 其 边界 是 C。 方程 (7.23) 表明 涡 管 的 截面 
RSA QI SHRRSET ARERR MHA, BER 
理 说 随 着 流体 运动 的 环 量 是 常数 。 所 以 得 计算 DPIVDr。 但 


Pr D $ u-d=¢ Mat fu: du, 
Dt Di 。 c c Dt C 


式 中 的 如 是 沿 围 线 的 速度 矢 的 微分 量 ， 不 过 u du-d (Huu); 


所 以 如 果 速 度 场 是 连续 的 话 ， 右 侧 第 二 个 积分 等 于 需 。 用 方程 
(7.9) 或 (7.16), 并 引用 斯 道 克 斯 定理 


pr __ | [v x (+ vp) -dÀ = 1 cr XVs)- dA (7.24) 


EE, R E ERU RE SE BE FSB, TS EAR BY 
任 一 闭路 的 环 量 必 是 常数 .这 是 开 耳 芬 (Kelvin) 定理 
不 过 在 一 般 的 场合 , 郎 使 是 理想 气体 , 环 量 出 不 是 常数 ， 这 和 
密度 不 变 的 不 可 压 流体 的 情况 是 大 不 相同 的 . 事实 上 ,方程 (7.24) 
可 以 写成 
DL | [ve x v+] aA, (7.25) 
4 p 


FB 7.25 TUTE DA, MOR 2 是 常数 ( 璧 如 不 可 压 流体 ), BRT — 
般 地 总 是 常数 ; 对 于 可 压缩 流 来 说 , 方程 (7.25) 可 以 作 如 下 的 解 
FE: 在 流 场 里 作 等 间距 的 p = const 曲面 族 和 1/o = const 的 曲面 
族 ， 这些 曲面 围 成 一 系列 的 管子 .党 着 随 流体 运动 的 围 线 C 的 环 
量 ,在 单位 时 间 内 的 变化 是 正比 于 C 所 包 的 管子 数 旦 的 ， 这 条 定 
理 是 伯 耶 克 纳 斯 (Bjerknes)n 提出 来 的 . 

鉴 点 量 . 在 本 节 讨 论 中 引进 r 为 单位 质量 的 给 和 动能 之 和 ， 
或 总 能 量 。 这 个 量 通常 称 为 流体 单位 质量 的 “ 驻 点 痊 值 ?, 其 概念 
是 这 样 的 : 用 某 种 办 法 从 流 场 任 何 一 点 取出 流体 ， 写 具有 在 这 一 - 
点 的 当地 压强 、 密 度 、 温度 和 流速 , 然后 以 减低 流速 的 办 法 把 这 团 
流体 作 等 炳 压缩 , 直到 流速 降 到 零 为 止 , 郎 达到 沸 止 为 止 ， 流 体 微 
团 在 滞 止 条 件 之 下 所 达 的 最 后 状态 ， 称 为 原 取 出 点 的 流体 的 驻 点 
状态 。 SOWA, BREE, PARMAR BESS 
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J. YARATAR, RELENST SE. WA SREB 
W, FE Ya BIE RAS RA TRS ED tt 4 BH, RAKES, BOA 
WHERE, EREE, SE, RE, MPS 
SCOG STIS AY Fe ee SEO Ge, BABE RAE Ls a 
量 。 在 亚 声速 流 中 ， 用 皮 托管 量 得 的 压强 只 是 驻 点 压强 的 近似 值 
ME. 在 超声 速 流 里 , 电 于 皮 托 管 前 存在 激 波 , 所 测 压强 较 驻 点 压 
RES, 

FECA HRM AK, TRESE ZAKRES 
X., OM EmA, EE SHR 4 RL, 

|u|? 

rp/p) ` 
BE (1.5), (5.7), (5.8), (5.9) 和 (7.11), MBE cp 是 常数 ， 
就 有 


M 一 


to mM, (7.26a) 


PR Je FAD BT HA SES at FA BS 5 a TA ASE, BAT 
(7.26b) 


在 完全 气体 , AAS AE BIEL, EOE A RS SEE 
A BAL, 

当地 马赫 数 很 低 的 时 候 , 方程 (7.26) 可 以 展 为 M 的 级 数 , 例 
如 


Logt mt mt CTN yoy ..., 

p 2 8 48 
所 以 

Opa 2 tye 427 pt... 

pP—p zelul (+++ 元 M 小 (7.27) 
驻 点 压强 和 静 压 之 差 p 一 p EJERE., 4MB hit, Hee 
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近似 地 等 于 + p|al:， 这 个 量 常 称 为 动 压 , 


A.8. 理想 气体 的 绝热 流动 . 伯 努 利 方程 


在 前 节 的 讨论 里 已经 说 过 ， 在 航空 方面 的 大 多 数 气 体 动 力学 
问题 , 除 边 界 层 外 , 粘性 和 热传导 的 作用 都 可 以 略 去 不 计 。 此外， 
除 有 关 燃 伐 的 问题 之 外 ， 和 外 加 热量 不 是 零 便 是 很 小 。 押 以 在 普通 
条 件 之 下 ,气体 很 象 理 想 气 体 ， 因 此 ,气体 动力 学 的 基本 问题 之 一 
是 研究 理想 气体 的 绝热 流动 ， 本 节 以 下 的 讨论 大 部 分 是 关于 这 个 
问题 的 . 

设 流 动 是 无 旋 的 ,因而 有 一 个 速度 势 P(xi, t) 存在 ， 

u = Vġ, (8.1) 
WMQ=Vvxu=V X (Vd) AHAB. 动力 学 方程 (7.13) EK 
a 1 d 
v (Fe + lup |£) -=0, 
这 就 是 说 括号 中 的 量 不 是 一 个 空间 两 数 , 但 可 能 是 时 间 的 画 数 ;这 
个 时 间 夯 数 可 以 包 藏 在 + 的 定义 里 , 而 不 影响 方程 (8.1) 所 示 关 
系 . 所 以 


So 十 二 [aa 十 | 2 = const, (8.2) 
Or 2 p 


这 是 非 定常 无 旋 流 的 伯 努 利 方 程 ， 可 以 看 作 是 动力 学 方程 的 第 一 
次 积分 ， 寂 给 出 压强 或 密度 和 速度 势 之 间 的 简单 关系 。， 式 中 积分 
项 必须 按 等 炳 的 压强 -密度 关系 计算 , 因为 前 面 说 过 , 无 旋 流 通常 
只 能 在 全 流 场 中 单位 质量 的 炳 为 常数 时 才 得 保持 . 在 完全 气体 ， 
方程 (8.2) 变 成 


Ob Ljuj 4 7 Pp 10 = const., (8.3) 
Or 2 7 一 1 2 Or 


所 以 无 旋 流 工 不 一 定 要 均 能 , 朗 全 流 场 各 点 的 总 能 如 不 一 定 相 
同 。 只 在 定常 流 , 朗 06/0: = 0 时 ,无 旋 流 才 是 均 能 的 . 

如 果 流 动 是 有 旋 的 话 , 也 就 是 涡 量 不 为 零 的 话 ,一 般 的 第 一 次 
积分 不 存在 ， 只 有 流动 定常 时 ， 第 一 次 积分 才 存 在 .如 前 所 示 在 
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绝热 流 中 , 沿 流体 微 团 的 迹 线 或 沿 流 线 的 总 能 如 是 常数 ， 但 不 同 
的 流 线 可 能 有 不 同 的 数值。 


eT P=h+ > |u|? = const., (8.4) 


这 实在 是 个 能 量 方程 ， 但 在 形式 上 和 定常 无 旋 流 的 伯 努 利 方程 一 
样 ,因为 后 者 是 均 能 的 ， 二 者 的 差别 当然 在 于 ,定常 的 均 能 流 不 一 
定 是 等 箭 的 ;而 定常 的 无 旋 流 都 必定 是 等 箭 的 。 举 个 例子 来 说 , 具 
有 膀 体 旁 曲 激 波 的 绕 物 体 的 定常 超声 速 流 是 均 能 的 但 不 均 炳 . 
以 上 的 讨论 表明 流动 有 四 种 特征 , 序 定 常 性 、 有 旋 性 、 均 业 性 
和 均 能 性 。 按 这 些 特征 来 划分 ， 理 想 气体 的 绝热 流动 有 如 下 所 示 
的 八 种 可 能 ,“ 十 ”号 表示 具有 表 内 所 示 特 征 ,“ 一 ”号 表示 没有 那 种 
特征 。 例 如 在 定常 性 下 的 “十 ”号 表示 流动 是 定常 运动 ,在 有 旋 性 
下 的 “一 "号 表示 流动 是 无 旋 运 动 . 
类 型 EEE AiE SME ” 均 能 性 
+ - + 


l 
1 + 上 十 


on nN hh WwW NH = 
1++1+1 
二 十 十 十 十 | 


一 + 


VAM Rite, 如 为 定常 的 则 属 第 一 类 ; 非 定常 的 属 第 2 类 . 
流 过 机 村 的 跨 声 速 流 ， 如 为 定常 的 则 属 第 3 类 ; 非 定常 的 属 第 4 
类 。 设 绕 物 体 定常 流 的 自由 流 是 非 均 匀 的 , 即 切 变 流 ,如 果 自 由 流 
的 温度 和 压强 是 均匀 的 但 速度 分 布 不 均匀 , 那 末 属 第 5 类 ;否则 属 
第 6 类 .如 果 此 种 流动 是 非 定常 的 , 则 属 第 7 或 第 8 类. 

如 果 运 动 的 方程 有 第 一 次 积分 存在 ， 那 未 流 场 的 数学 分 析 的 
复杂 性 大 为 减少 ， 上 面 说 过 ,无 旋 流 或 定常 流 就 存在 第 一 次 积分 . 
YARRA EAC, AM eee. 将 于 以 下 
各 节 详 细 讨 论 这 几 种 流动 。 
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A.9. Cit. BRS 


SUSE FOREN , VA TRE BH p(xi, 1) FE, MA BID 
PBB BA ABA (8.2), ARERR KE ER BT 
以 能 量 方程 和 状态 方程 就 给 出 三 变数 p, p, T 之 中 任意 两 个 之 间 
BUSES. 从 等 炳 关系 和 伯 努 利 方 程 可 以 得 到 一 个 有 速度 势 和 
包括 热力 学 变数 p, p, 了 之 一 在 内 的 方程 。 例如 ， 比 热 为 常数 的 
完全 气体 ,如 果 远 离 物体 处 的 流动 是 均 与 定常 流 , 且 远离 物体 处 的 
驻 点 条 件 也 是 均匀 的 , 合 驻 点 变数 注 以 上 标 0, 则 


rY-1 
-yT 0 
2 44 uj +—t_ £(2) ”一 一 了 £ 
Oz 2 7—1po\p y—i p 
= @Y, (9.1) 
9G 4 Ljuj + 7 r(e = # 
Oz 1 p°\ 9° 7 一 1 æ 


=— (wF, (9.2) 
7 一 


be 4 Lyyp+—_#(Z)-_2 
Or 2 y¥—1 o o 


一 一 上 (æF, (9.3) 
-一 


Ob 1 2 1 2 = 1 O42 
SP 4 + = 9.4 
a tg ht ee, (94) 


式 中 的 “ 是 声速 , 定义 见方 程 (1.4)。 当 速度 势 , 从 而 lul 都 确定 
之 后 ,利用 上 式 就 可 以 计算 压强 、 密 度 和 温度 . 
于 是 问题 归结 为 求 速度 势 . 因为 伯 努 利 方 程 是 积分 动力 学 方 

程 和 能 量 方程 的 结果 ， 所 以 唯一 没有 用 过 的 运动 方程 是 连续 方程 
(7.7), 该 方程 可 以 写 为 

43e ry.ut 工 a.vp 一 0. 

p Of p 
{AEN EM, AR Ae e ER, 上 式 也 可 写 为 
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(+22 + u-lyp)tv-u=0, (9.5) 


elie e Of 
而 一 般 情况 下 , dp/dp 只 是 和 声速 有 关系 , 即 方程 (1.3) 
Mp 
dp 
将 伯 努 利 方程 对 时 间 : 微分 ， 
Po yu. 2u, (9.6) 
ðr? Or P Of 


从 这 个 方程 可 以 解 出 (1/p)(6p/6:)。 将 此 结果 代入 (9.5) 式 , H 
用 动力 学 方程 消去 (1/p) Ve, 方程 (9.5) 最 后 变 成 
_10¢ u ôu 


Vv: 一 一 一 十 2 一 tae u-Vo)u 9.7 
wate a [Cu.v)m (9.7) 
或 

1 Op 2 Od 1 Od 
Vo 一 一 十 二 十 — u;uj > 9.8 
PT | a Don a” x,0x,; (9.8) 


到 可 以 用 以 由 表示 的 伯 努 利 方 程 来 计算 . 对 于 完全 气体 来 说 ,oa 
由 方程 (9.4) 给 定 。 如 此 说 来 方程 (9.8) 可 以 看 作 是 仅仅 一 个 未 知 
量 4$ 的 方程 .这 较 之 原始 运动 方程 组 是 一 个 大 简化 . 
当 声 速 和 流体 的 速度 比 起 来 十 分 大 的 时 候 , 则 方程 (9.8) 的 右 
MWh Oe, ORAS, FE 
Vp =0, a>, (9.9) 
这 是 不 可 压 流 速度 势 的 微分 方程 ， 从 数学 上 说 来 , 方程 (9.9) 和 一 
般 方程 (9.8) 根 本 不 同 , 方程 (9.9) 是 未 知 落 数 > 的 线性 方程 , 而 方 
程 (9.8) 却 是 非 线性 的 。 所 以 解 方 程 (9.8) 的 数学 问题 远 较 解 相应 
的 不 可 压 流 问题 困难 得 多 。 从 物理 上 看 , 方程 (9.9) 表 明 无 旋 的 不 
可 压 流 的 速度 分 布 问题 可 以 不 用 伯 努 利 方 程 求解 。 伯 努 利 方程 是 
在 速度 决定 之 后 用 来 求 压强 p 的 .事实 上 , 方程 (9.9) 可 以 从 连续 
方程 和 无 旋 条 件 去 推 得 . 所 以 在 无 旋 的 不 可 压 流 ， 速 度 场 仅仅 为 
运动 学 条 件 所 决定 ， 和 动力 学 问题 完全 没有 关系 ， 可 压缩 流 就 不 
是 这 样 。 此 处 运动 学 问题 和 动力 学 问题 紧密 地 连结 在 一 起 ， 因 为 
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速度 势 的 方程 , 朗 方 程 (9.8 ), 必须 与 伯 努 利 方 程 同时 求解. 

可 压 流 和 不 可 压 流 的 另 一 重要 区 别 是 所 谓 的 “ 表 观 质量 “ 问 
题 ， 这 是 物体 在 远 处 静止 的 流体 中 作 加 速 运动 时 表现 出 来 的 一 个 
量 .在 不 可 压 流 ,任何 一 瞬时 的 流 场 只 为 运动 学 的 条 件 所 决定 , 故 
物体 加 速 运动 所 引起 的 瞬时 流 场 和 以 此 瞬时 速度 作 定 常 运动 的 流 
场 完 全 相同 。 只 有 在 计算 物体 表面 压强 时 ， 运 动 的 非 定 常 效应 才 
通过 伯 努 利 方程 里 的 06/0: 项 显示 出 来 ， 但 中 是 和 瞬时 速度 成 
正比 的 ;所 以 加 速度 的 附加 压强 效应 与 速度 的 时 间 变 化 率 或 物体 
的 加 速度 成 正比 。 积分 物体 表面 压力 即 得 作用 在 物体 上 的 合力 . 
所 以 由 于 加 速度 而 来 的 附加 力 是 正比 于 物体 的 瞬时 加 速度 的 。 结 
果 使 物体 四 周 的 流体 加 速 的 效应 可 以 看 作 是 有 一 额外 质量 和 物体 
一 起 运动 。 此 附加 质量 常 称 之 为 表 观 质量 ， 宪 是 一 个 和 加 速度 无 
关 的 常数 ， 而 只 决定 于 物体 的 形状 、 这 就 大 大 简化 物体 在 流体 中 
作 加 速 运 动 的 计算 .从 物理 上 看 来 ， 这 个 简化 实际 是 扰动 在 不 可 
压 流 体 里 瞬时 传 肖 全 场 的 结果 .因此 加 速度 的 效应 完全 为 该 物体 
的 瞬时 加 速度 所 决定 ， 在 可 压 流 中 ,扰动 的 传播 速度 是 有 限 的 ; 运 
动 的 历史 工 不 能 立 部 不 起 作用 .因此 在 某 一 瞬时 ， 加 速度 的 效应 
不 仅 为 该 瞬时 情况 所 决定 ， 而 且 为 前 此 一 柔 烈 的 情况 所 决定 。 于 
是 加 速度 的 效应 就 和 物体 的 瞬时 加 速度 很 少 有 关系 了 ， 而 表 观 质 
量 的 概念 就 没有 什么 用 处 . 

如 果 运 动 是 定常 的 ， 所 有 对 时 间 的 导数 都 等 于 零 . 这 时 速度 
势 的 方程 (方程 9.8) 变 成 


2, 1. yp 
Vd = a @x,0x;° (9.10) 


声速 a 仍 可 用 伯 努 利 方 程 计算 ， 现 在 该 方程 中 去 掉 了 06/0: W, 
对 二 维 流 可 以 分 x 二 xj, y = x; u = th, v = a, 方程 (8.10) 变 成 


u’\ Od uv Op v\ Od _ 
(1-4) 58— 2 oe + (1-8) S$ o) 


*“ 表 观 质 量 ” 在 许多 苏联 书 上 称 为 附加 质量 一 一 译 者 注 . 
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A.10. 无 旋 流 的 变 分 法 


无 旋 流 问题 除了 从 微分 方程 (9.8) 求 解 速度 势 之 外 ,还 有 可 
能 形成 一 个 变 分 问题 ， 所 谓 变 分 问题 是 由 的 基 个 积分 在 中 作 微小 
Ay de aE te (Virtual variation ) 时 应 该 有 个 逗留 值 (Stationary value), 
这 样 一 种 作法 首先 由 贝 特 门 (Bateman) 提出 ， 后 来 为 王 启 德 
(C. T. Wang)" 所 推广 , 

为 了 明确 起 见 ， 考 虑 求解 包围 在 固定 固体 空间 边界 中 的 一 种 
JEE MACK, FERN Wt =n, 在 这 个 空间 边界 内 的 
速度 分 布 是 给 定 的 . FAT NW: = ,空间 边界 内 的 速度 分 布 
也 是 给 定 的 . Ea o FERRERIAS E, 而 是 必须 
HEA, BERR i A NARRATE ST, E t Mt i AR 
度 就 是 给 定 了 , 那 末 在 所 和 妃 时 边界 以 内 空间 每 一 点 上 的 中 本 身 
的 值 除了 一 个 常数 之 外 也 是 给 定 的 。 可 以 取 常 数 为 零 ， 而 不 致 改 
变 问 题 ， 换 旬 话 说 , 问题 已 变换 为 这 样 : E n M n ito ERRET 
的 ,同时 (64/6xi) ni 在 边界 上 必须 永 为 需 ， n 是 边界 的 单位 法 向 
RK. 现在 考虑 积分 


I= (if pdx, dx dxs di, (10.1) 


HRASRTRALAANAMFAMA =n RSD t—n, PÈ 
压强 ,从 伯 努 利 方 程 (8.2) EERS STH, WINK, tit 
动 问题 等 于 求 这 个 积分 了 的 逗留 值 ,部 

6I = 0, (10.2) 
因为 伯 努 利 方 程 是 从 动力 学 方程 用 了 速度 势 $ 推导 的 结果 ， 无 旋 
条 件 已 经 是 满足 了 的 ， 唯 一 的 没有 用 到 的 方程 是 连续 方程 (2.5). 
要 证 明 贝 特 门 的 说 法 就 是 要 证 明 方 程 (10.2) 和 连续 方程 (2.5) 是 等 
价 的 。 利 用 伯 努 利 方程 (9.1 ), 方 程 (10.1) 成 为 


i 


X dx dx dx3 dt, (10.3) 
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ee 


MZ 6 — SEES Od, M 


n- -ejb -e C9) 


x (2 + ob 228) dx, dx, dx; dt. 
Oz Ox; Ox; 


按 (9.2) 形 式 的 伯 努 利 方 程 看 来 , 这 个 积分 的 被 积 式 的 第 一 个 因子 
等 于 密度 p。 故 


sf = 一 MIE (es + e £32) dr, dx; dx; dt, (10.4) 
Ox; Ox; 


9 Ob 288 9 (0% 98) — s 2- (p $) 
P Ox, Ox; Ox; $ $ Ox; P Br 全 


所 以 方程 (10.4) 可 以 作 如 下 变换 ,第 一 项 对 # 作 分 部 积分 , 第 二 项 
用 方程 (3.1) 的 斯 道 克 斯 定理 作 变 换 。 结 果 是 


ôI 一 一 {f (dp 28 dx, dxa dx3 一 {a | dA(5Puini) 


+ fas( 十 Zeti) dy dr des dt, (10.5) 


式 中 的 面积 分 是 沿 固体 边界 表面 进行 的 。 第 一 个 积分 为 零 ， 因 为 
aM t EFAG p eae TH, 所 以 在 二 和 刀 时 ,一 切 空间 点 的 54 者 
AS, 第 二 个 积分 也 是 零 , 因为 在 边界 的 表面 上 umn 一 0。 所 以 
方程 (10.2) 的 定义 不 过 是 方程 (10.5) 的 第 三 个 积分 必须 等 于 雾 而 
已 。 但 除了 在 志和 如 时 之 外 ,59 是 任意 的 ; 故 方程 (10.2) 只 有 在 
连续 方程 得 到 满足 时 , 朗 

Bo 4 domy 

Oz Ox; 
时 , 才 可 能 成 立 。 这 就 证 明了 贝 特 门 的 说 法 . 

至 于 定常 流 问题 , PAVE t Ml NY 由 ,积分 工 简化 成 为 


了 一 {ff p dx, dx, ax3, (10.6) 
式 中 2 是 用 伯 努 利 方程 作 计 算 的 ,当然 该 式 中 的 68/61 项 是 不 存 
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ERS, 

变 分 方程 (10.27 的 应 用 就 是 用 直接 的 瑞 利 〈Rayleigh) 与 李 兹 
(Ritz) 的 变 分 法 去 计算 流动 。 要 点 是 选取 一 个 画 数 近 似 地 表达 呐 
正 的 速度 势 p, 所 选 画 数 要 满足 边界 条 件 以 及 上 一 志和 上 一 所 时 
的 条 件 ， 而 且 还 得 包含 几 个 任意 常数 .运算 的 过 程 包含 把 近似 丙 
数 代入 方程 (10.1) 或 (10.6) 里 去 , 视 非 定常 流 还 是 定常 流 而 定 , 然 
后 积分 之 。 这 样 一 来 ,积分 了 就 成 了 近似 画 数 的 未 定常 数 的 画 数 . 
最 后 用 方程 (10.2) 定 这 些 常 数 . 这 样 定 得 的 含有 常数 的 丽 数 就 是 
里 实 速 度 势 的 近似 表达 式 . 

王 启 德 货 发 现 应 用 变 分 法 于 无 限 空间 的 定常 流 必 须 克 服 两 重 
困难 ， 第 一 个 困难 是 ,这 时 积分 (10.6) 是 无 限 大 ,因而 上 述 变 分 问 
题 的 处 理 是 不 恰当 的 。 第 二 个 困难 是 , 方程 (10.5) 里 的 面积 分 现 
在 不 等 于 替 ， 因 为 曲面 本 身 是 无 限 大 。 第 一 个 困难 在 亚 声 速 流 这 
样 克服 : 在 方程 (10.1) 中 彼 积 式 不 用 了 而 用 ? 和 对 应 的 不 可 压 流 
中 的 Pine 二 者 之 差 , 即 用 


r = |ÈÈ Co — Pine) dn dr dss, (10.7) 


如 果 不 可 压 流 的 解 是 已 知 的 , 那 末 Pine 的 积分 是 个 固定 的 值 :所 以 
用 差 来 代 p 并 不 改变 上 面 关于 OF 的 推理 . 但 这 样 的 技巧 就 使 那 
个 积分 变 成 有 限 的 了 .第 二 个 困难 这 样 克 服 ， 用 速度 势 的 近似 画 
数 去 实际 算出 (10.5) 式 中 的 面积 分 来 ,然后 用 它 去 修改 方程 (10.7) 
的 积分 , 使 改变 后 的 积分 之 变 分 恰 能 对 消 方 程 10.5) 中 的 麻烦 的 
面积 分 (详情 请 参看 [12]). 


A.ll. 绍 热 定常 流 


前 一 节 里 证 明 过 ,在 理想 气体 的 绝热 定常 流 里 , 沿 流 线 总 能 或 
EAGER, BK VY 必 是 一 个 与 流动 方向 垂直 的 向 量 . 但 
按 方 程 (7. 18) 看 ， 定常 流 必 有 

u xX Q = FR — TVs, (11.1) 
PORE (VA — TVs) 是 垂直 于 流动 方向 的 。 不 过 , V 如 已 经 垂 
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ETR TS PA Vs Be PA i, SA 
是 常数 . 激 波 之 类 的 突 跃 变化 当然 是 例外 这 种 突 跃 面 出 现时 ， 
灶 也 作 突 跃 地 变化 。 因 此 ,在 激 波 与 激 波 之 间 , 沿 一 条 流 线 的 定常 
Sy PCIE SN, BRS, 这 时 流 绪 
BLA BSE Rh p FUME sR FASE, TF “Re AEE” 
Umax IRPJ, Uma 的 定义 是 


如 一 > tax, (11.2) 


所 以 um: ECR eB Uh BE SET BY A RHE SE, HEE a 
部 转变 成 动能 时 的 流速 。 当 然 沿 任何 一 条 流 线 wmax 都 是 常数 . 

现在 用 最 大 速度 um: Ru, 除 出 来 的 无 量 纲 向 量 记 为 W; W 
可 以 叫做 减 缩 速度 矢 : 


w=. (11.3) 
用 能 量 方程 (方程 8.4) 得 

Ži -pm 

pi Ww, (11.4) 


式 中 的 W? = W . W, 或 W 的 绝对 值 的 平方 。 设 om 为 驻 点 密度 
(其 沿 流 线 之 值 也 为 常数 ), 则 连续 性 方程 是 


o=: (on) = y [ee /2H (5) Ww) 
p 
= 0° 27 « (5 w] + (5) W: WPN ZR ), 


P W2 如 沿 流 线 是 常数 ; 所 以 写 的 梯度 必 垂 直 于 速 矢 或 W， 因 此 
上 式 右 侧 第 二 项 是 雾 ， 从 而 


p 
v | 2) w] = 0, (11.5) 
对 于 完全 气体 来 说 ,利用 上 方程 (11.4) 可 将 方程 (11.5) 写 为 
v- [Cl 一 也 W] 一 0. (11.6) 
现在 动力 学 方程 变 成 
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(W-Y)W+ 1 (1—-Ww)Vinp=0, (11.7) 
Y 


如 果 不 用 向 量 W ， 而 用 当地 声速 除 速度 矢 u, 或 对 完全 气体 
来 说 用 WYp/p REEK, 那 就 得 另 一 个 无 量 纲 向 量 , 序 马赫 数 向 
aM, 用 新 向 量 表示 , 则 完全 气体 的 连续 方程 和 动力 学 方程 成 为 


Y+ 
g- l(a +2 = l ut) =? M| = 0, (11.8) 


(M: vV)M— ZZL M(v-M)+tvinp=0, (119) 
ytl Y 


向 量 W 和 M Sake (Munk) 和 浦 黎 姆 (Prim), Foe 
(Hicks)、 顾 恩 特 (Guenther) FIFE (Wasserman)? 等 所 引用 
的 。 在 每 一 情况 下 ,连续 方程 和 动力 学 方程 总 是 解 两 个 未 知 数 W 
和 了 ?抑或 M 和 ?的 方程 组 .不 过 应 注意 ,在 方程 (11.6) 和 《11.7) 
或 (11.8) 和 (11.9) 中 驻 点 烩 都 没有 显 式 地 出 现 。 而 方程 (11.6) 
和 (11.7) 或 方程 (11.8) 和 (11.9) 的 一 个 解 确实 是 包含 了 相当 于 给 
Ay SRR AR BE SB Te, 但 每 一 族 解 的 ?是 一 样 
的 ,和 所 给 的 驻 点 烩 值 无 关 。 警 如 , 若 已 知 绕 某 给 定 物体 和 给 定 流 
5 TB) KD 7 YE FS De , PRAE EARED AAI, RR 
EDRR BOT I, NET PE RIZE 
AR, 同一 个 压强 场 仍 能 保持 流动 平衡 . 

对 于 完全 气体 ,方程 (11.1) 用 W 表达 则 为 

Wx(vyxW)_ 7 一 1 
1 一 w? 2y 
式 中 的 是 驻 点 压强 ,党 一 条 流 线 仍 是 常数 .就 方程 (11.10) 取 旋 
度 , 印 得 克 罗 柯 《Crocco) 方程 : 
W x (V X w] = 0 

1— W? ° 
KT AEF LEECH 9 Se AISA ET SE SY, 

作用 在 等 速 运动 的 物体 上 的 力 . 现在 应 用 第 4 节 中 的 积分 
定理 来 计算 一 下 作 等 速 运 动 的 物体 所 受 的 力 、 力 徐 和 热传导 ， 积 
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Vln p, (11.10) 


vx | (11.11) 


分 沿 包围 物体 的 控制 面 41 进行 。 计算 用 的 方程 是 方程 (4.10)， 
《4.11) 和 (4.12)。 如 将 控制 面 41 取 成 一 个 远离 物体 的 大 曲面 , 远 
到 使 压强 、 密 度 、 温 度 和 速度 在 那里 随 空 间 坐 标的 变化 已 十 分 组 
慢 ,这 些 量 的 梯度 是 如 此 之 小 以 至 于 在 41 上 的 粘性 应 力 和 热 通 量 
都 可 以 略 去 不 计 ， 这 也 就 是 说 , 在 4 面 上 , 边界 层 和 激 波 的 效应 
已 经 扩散 到 只 剩 下 物体 后 面 的 一 条 宽 宽 的 尾 迹 了 . 这 时 方程 
(4.10),《4.11) 和 (4.12 ) 变 成 


f 
i= \, [ pn; + pu;Cujnj)| dA, (11.12) 
1 
M; = |, Sijl pxing + pxiur( un1))] dA, (11.13) 
1 
H=— | o( + L um) un; dA, (11.14) 
dı 2 


式 中 的 F;, Mi 妃 分 别 是 作用 在 物体 上 的 各 分 力 、 绕 各 坐标 轴 的 
力矩 以 及 传 给 物体 的 热量 . 另外 还 有 一 个 从 方程 (4.2) 来 的 连续 
方程 
| pum: dA = 0. (11.15) 
现在 分 Pos Pos Tos ho 代表 物体 远 前 方 的 均一 压强 、 密 度 、 温 
度 和 单位 质量 的 烩 值 ， 也 就 是 未 经 扰动 的 各 量 。 合 忆 代 表 未 扰 自 
由 流 的 速度 (xi 方向 ), 并 记 上 述 各 物理 量 与 均一 值 之 着 为 5p,6ui， 
Sh 等 。 则 物体 的 阻力 或 Fi 便 是 


Fi 一 | (Po + 6p)n dA 一 | CU + ui)puini dA, 
1 Ay 


pom 在 闭 曲 面 41 LMR AR, 再 根据 连续 性 方程 (11.5) HU 
为 常数 的 事实 , U(puini) 在 曲面 41 上 的 积分 也 是 雾 ， 此 外 ;在 A 
上 ,流动 情况 和 均一 情况 相差 极 微 , 只 需 考 虑 一 次 微量 就 行 了 ， 结 
果 阻 力 Fi 成 为 
Fi 一 一 |, (8p + pU sum dA, (11.16) 
方程 (11.14) 经 类 似 的 变化 成 为 
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H= — | PoU(6h + Usum dA, (11.17) 
Ay 
从 方程 (11.16) 和 (11.17 ) 中 消去 sw 得 


0P 一 已 一 | paU (on — 22) m dA. (11.18) 
Ar Pœ 


如 so REKETA II—I eT, Os 代表 在 4, 曲面 
E s 和 so 的 差 , 那 末 现在 按 小 量 偏离 的 近似 计算 得 


Ta6s = 6h— 2 (11.19) 
代 方程 (11.19) 入 方程 (1.18)， 
UFI—H=T, | paUdsnid A. 
另 一 方面 ， 用 连续 性 方程 很 容易 证 明 下 列 的 关系 : 
| ps(ujn;) dA 一 L PoU ssm dA, 


故 最 后 得 
UFı— H = To D ps(un;) dA, (11.20) 


BRCAL20) SF AKRA (Oswatitsch)™ Hi, BER 
很 有 用 的 定理 : 阻力 与 物体 运动 速度 的 乘积 减 去 物体 所 吸收 的 热 
量 等 于 自 包 围 物体 的 一 个 大 控制 面 流 出 的 炳 通 量 乘 以 自由 流 的 温 
EB, 不 过 , 要 这 条 定理 成 立 ,控制 面 必须 取得 足够 大 , 大 到 能 把 所 
有 能 量 耗 散 过 程 尽 包 在 里 面 ， 如 边界 层 尾 迹 和 激 波 。 其 实 这 条 定 
理 是 按 流 动情 况 和 均匀 未 扰 流 相差 极 微 的 假设 推导 出 来 的 ， 故 不 
HE FAT RA ARSE Sep eit, FREAK EEN, RE 
记 这 一 限制 . 

在 无 粘 无 热传导 、 无 激 波 的 理想 可 压 流 里 , 物体 所 吸收 的 热 
BAST, KA RWMRE SN RRA, BOK 
种 场合 , RHE (11.19) 的 近似 程度 说 来 ,阻力 是 没有 的 ， 这 就 是 
所 谓 的 达 明 培 疑 题 ， 是 西 负 道生 (Theodorsen )ng 首先 提出 讨论 可 
压缩 流 的 达 朗 培 疑 题 的 。 不过， 在 跨 声 速 流 和 起 声速 流 里 总 有 激 
波 出 现 的 ,所 以 对 可 压 流 来 说 , 达 朗 培 疑 题 远 不 如 在 不 可 压 流 里 那 
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RAD, BROT ACA ERRES , Plane LRN = See Bh R 
Uh SN FAC 11.19) HRA ST , BD AR, = HL 
BA Dt Re 


A.12. 二 维 定常 均 能 流 ， 流 画 数 
Mt HR th eR, 流动 可 以 是 定常 均 能 的 , (IS, BA 
如 绕 跨 嘻 速 飞行 物体 并 具有 离 体 激 波 的 流动 就 是 如 此 ， 这 时 全 流 
ASE RIAA EP BR, SEF eee. (Bie 
BME aS, FAD PEC11.16) (11.11), 4 


0=V x {ra — WHIWIx [| 


| vxW 


es Y| a- OW 


7 一 vx W 
TIG OWN ap 


— pyar) .vxW 

十 《1 一 Ww?) w |v aa (12.1) 
二 维 流 的 Y x W 只 有 一 个 分 量 , ERE TR, PAE 
AWABBIMPRRETE, BA, WR « 和 wv 是 W 在 x 和 
y 方向 的 两 个 分 量 , 则 方程 (12.1 ) 变 成 

Ov Bu 

u 4 yp 8 Ox oy dy | 

( ðs a) la 一 al °. (12.2) 
PRE SS FB iit BO tee 4a R BF a W 和 物理 流速 向 量 u 的 差异 在 于 一 
个 常数 因子 . 所 以 方程 (12.2) 可 以 解释 为 涡 强 除 以 (1 一 w/o) 
所 得 的 商 沿 流 线 为 常数 .只 要 没有 激 波 之 类 的 突 跃 面 存在 , (1 一 
Ww?) RAY HEE RR IEC, PA, 在 定常 二 维 均 能 流 中 ,党 任 
何 一 根 流 线 ， 涡 强 是 直接 与 当地 压强 成 正比 的 。 这 条 定理 是 克 罗 
AYE HSE BOM), 

按 克 罗 柯 的 办 法 定义 流 画 数 为 
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Ob ull — WM), (12.3) 
Oy 


Ob = — y — 21y —1) 

Bw (1 一 w’) . 
这 样 一 来 连续 性 方程 (11.6) 自 动 满足 ” 流 画 数 为 常数 的 曲线 是 流 
Ay dn 代表 两 条 相 邻 的 流 线 几 和 由 十 dy 之 间 的 距离, 则 由 方 
程 (12.3) 得 


db = (1 — WYO W dn = £ W dn, 
p 


因此 ，pW dn 或 ou dy 是 和 流 线 之 间 的 质量 流量 成 正比 的 。 但 各 
SPEER s 不 同 , 不 同 流 线 上 的 o 也 各 不 相同 。 所 以 dy 本身 
不 直接 代表 上 几 和 4 十 44 两 条 流 线 之 间 的 质量 流量 的 。 这 样 一 
来 , 当 流 动 越过 激 波 时 ,由 于 炳 和 驻 点 密度 p? 有 突变 , ERY EE 
有 突 跃 。 方程 (12.3) 所 定义 的 流 画 数 有 突 跃 就 可 能 造成 不 便 。 陶 
ARBACTollmien )°7 试图 补救 此 缺点 ,另外 引用 一 个 流 汞 数 4 ,4' 和 
克 罗 柯 的 有 如 下 的 关系 
dp’ = odd, 
SRK, y 在 越过 激 波 时 就 连续 了 , 但 J' 的 方程 却 十 分 复杂 . 
前 一 段 已 经 说 过 ,方程 (12.2) 意 味 着 
Ov Ou 


Ox Oy 
GQ — WIAD 
WY PARE f(b) 的 物理 意义 可 以 就 现在 这 个 问题 把 方程 (11.1) 
分 析 一 下 。 设 7 是 流 线 的 法 线 ; 当 z 与 速度 矢 的 相对 位 置 和 ? 与 
x 的 相对 位 置 一 样 时 , 则 zz 是 正 向 。 那 未 方程 (11.1) 给 出 


= fy). (12.4) 


而 按 流 画 数 的 定义 (方程 12.3) 


wi — WD = 2e | 
dn 


e 47 a 


Ov Ou 
Ox Oy — _ 1 ds 
G — wan b= Fo (12.5) 


所 以 f(b) MARIERROTA., WE Nii ES 
的 , 则 ds/ay 为 需 , 这 时 涡 必 为 需 , 流 动 是 无 旋 的 . 

用 流 画 数 代 禁 方程 (12.5) 中 的 分 速 导 数 ,得 由 的 方程 如 下 ( 详 
情 参 看 [91): 


2) Op uv Od ( “) Op 
1—-4)SP 7 oe 十 (1 一 二) 一 
=) 0x’ ce AxOy ec/ Oy 


(1 ror- ( W? _ 
= a= wyrwa (Hie), 026) 
比 ,或 


c? = a — y 一 1 (1 _ w?), (12.7) 
Uraax 2 
RUF FEC 12.345 
xi — wrawr-» — (3y DY 
wi — w?) (24) +(#). (12.8) 


当 给 定 炳 的 梯度 ds/ay 之 后 ,譬如 说 由 曲线 激 波 给 定 , 方程 (12.6) 
便 是 瀛 画 数 的 一 个 非 线性 方程 . 

变 分 法 。 和 第 10 节 讲 的 无 旋 流 问题 一 样 ， 上 面 讨论 的 二 维 
有 旋 流 也 可 以 用 变 分 法 求解 ， 拿 下 烈 的 积分 来 看 : 


I= 人 eWay 一 zoom 人 1 十 工 一 w?) dxdy, (12.9) 
7 一 
A 


式 中 的 W? 由 方程 (12.8) 给 出 ,积分 域 4 为 流 线 及 与 流 线 相 垂直 的 
曲线 所 围 成 ， 变 分 法 的 基础 是 使 5 在 流 画 数 少 作 一 切 虚 变化 时 
MAR, Ep 

ôI = 0, (12.10) 
这 和 使 上 满足 微分 方程 (12.6) 或 更 原始 的 方程 (12.5) 是 完全 一 样 
的 。 这 个 说 法 的 证 明和 如下: 
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af 


ôI = {| er 一 一 a-wy™ (1 一 ztl w?) ow? 


ds 
—(1- p? “(i + cttw) t og| dx dy 
(1— w>’) > 二 1 R dp ° 
(12.11) 


但 流 画 数 的 定义 (方程 12.3) 给 出 下 烈 的 关系 式 


ap) _ ul — Wve) 一 一 (1 一 WYO sy?, 
Oy 7—1 


ZOL) gy — wave — 2 (1 — waver a? 
Ox y—i 


第 一 式 乘 以 u, 第 二 式 乘 以 v, 然后 相 加 部 得 
VELO typ) 5 (080 _ ,059 
ow Ty (1 Zt) 2(u By se), 
所 以 方程 (12.11) 可 以 写成 
DT = e™s OVA 2r EL —y Osh 
Jf om [22 ( ) 


7 一 1 Oy Ox 
— — g2) y+1,.\_1 ds 
(1 =w?) (rH r)i dh 54| ax dy, 
这 个 方程 可 以 用 分 部 积分 变 成 


GT 一 ff e7 OVA | 一 一 (24 — %) 
F y— 1i \Ox Oy 
EE (, a6 _ , 08) 

7 一 1 Z db Oy Ox 


1 ds 2 von ( ytl | 
一 一 -一 一 (LI -W 1 十 —-——- W Oldxd 
Z ap. ) 71 pdxdy 


— eye (_ 27 
$, =- — —1_) aW al, 
线 积分 是 沿 4 域 的 边界 线 进行 的 ， 而 这 个 边界 线 是 由 流 线 及 与 流 
线 相 垂直 的 曲线 围 成 的 ， 在 流 线 上 , Ob AB, 在 与 流 线 相 垂直 的 
曲线 上 ，W … 4dL 等 于 零 . 所 以 围 道 积 分 为 需 。 面积 分 中 的 导数 
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Of/Ox, Ob/Oy 都 可 以 通过 方程 (12.3) 表 为 速度 .于 是 


sr = 2I | RRON |2- — Ou. 
y—14 Ox Oy 


— E (a WD 2] dbdxdy, (12.12) 
2cp dd 


因 sp 在 A REER, MAF RSET S, 所 以 
(12.10) 式 条 件 要 求 方程 (12.12) HES PORTER AW, 
可 见方 程 (12.10) 是 与 方程 (12.5) 和 (12.6) 等 价 的 ， 

最 近 ， 克 有 罗 柯 吧 为 二 维 定常 不 均 灶 然 而 均 能 的 流动 引进 一 种 
PNB MR p。 9 的 定义 如 下 : 


2 2 2 
po 一 一 9， pu = OE, p +o = È, (12.13) 
y xoy 


把 这 些 关 系 代入 动力 学 方程 ， 动 力学 方程 恒 等 地 满足 。 9 的 好 处 
在 于 它 越 过 激 波 时 是 连续 的 。 不 过 从 连续 性 方程 推导 出 来 的 9 的 
微分 方程 却 是 个 非 线性 的 三 阶 方程 . 

轴 对 称 流 。 对 于 轴 对 称 流 ,最 方便 的 坐标 柔 是 * 和 7 的 柱 坐 
标 系 (x 沿 对 称 轴 , 7 与 此 轴 垂 站 )。 命 减 缩 速 度 W 在 x 和 7 两 个 
方向 的 分 量 为 x 和 wv, 那 末 流落 数 可 以 这 样 定义 (参看 方程 12.3) 


rull 一 WYD = -00 roll — Wye» = Ob (12.14) 
Or Ox 


克 罗 柯 发 现 ,现在 在 两 突 跃 面 之 间 , 沿 一 条 流 线 , 涡 强 正比 于 当地 
压强 * 及 7 的 乘积 .现在 涡 的 方程 是 


Oe _ du 
Ox Oy __ 
Fa warns T fe). (12.15) 
流 画 数 的 方程 是 
1 — EE) 2 Od A Pb 1 Ob 
c?) Ox? “ce Ox, 一 1-5 B77 r Or 
7 


= pa WIM (H Ko, (12.16) 
RRRS Cb) 和 梯度 ds/ay 的 关系 仍 是 方程 (12.5)。c RHR 
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(12.7), SREP A EEC 12.8 Mae EF SAMAR 
_ ry A q, (Ob) 
r WC 一 WOD (2) + (2). (12.17) 
本 问题 也 可 以 用 下 烈 的 积分 形成 一 个 变 分 问题 


I= (f eTA (a — WAY (1 + rt- w?) dx dr, (12.18) 
7 一 


式 中 的 WW 和 几 的 关系 见方 程 (12.17 ) ,积分 域 为 4，4 域 由 流 线 及 
与 流 线 相 垂直 的 曲线 所 围 成 。 方程 (12.15) 或 (12.16) 等 于 要 求 积 
Sy I FER > (ERIE AE Bh. 

非 均 能 流 . bm AE KEBAB ERODE ATR AD 
地 推广 到 非 均 能 流 . EXE WR ee W 及 其 分 量 u,v 的 
话 , 那 末 流 画 数 的 定义 方程 (12.3), 上 的 方程 (12.6) 以 及 辅助 方程 
(12.1) 和 (12.8) DAZ. ERZ Kb) 就 不 是 如 方程 (12.5) 所 
示 , 替 代 写 的 是 由 方程 (11.10) 而 来 ,现在 变 成 

—1dinop° 
Kg) = -L a (12.19) 

轴 对 称 流 也 是 这 样 . 当然 在 非 均 能 流 里 ， 不 同 流 线 的 最 大 速度 
umax 是 各 不 相同 的 ， 而 减 缩 速 度 W REDE u 和 vw 与 物理 流速 
下 及 其 分 量 的 差别 也 不 是 一 个 常数 因子 了 。 如 此 说 来 , 不 能 将 
《ax/azxr) — (8x/8y) 理解 为 与 流体 的 涡 量 成 正比 的 量 了 ， 


A.13. 二 维 定常 无 旋 流 . BEMER 


WORE ds/dy AS, 流动 就 是 无 旋 的 均 炳 的 。 这 时 方 
程 (12.6) 的 右 侧 是 零 ， 这 时 由 于 方程 (11.4) 的 关系 而 得 


(1 — wave 一 S (13.1) 


式 中 的 0° AGE REE BE, 0° SER MIG ARE HK, 而 且 在 整 
PRD BE HR, FFA FEC 13.1), HER BAS ELE RR 


Op d Op p 
R Lu, ZE, 13.2 
By P” Ox P” (13.2) 


这 时 ， 两 条 流 线 的 罗 值 之 差 就 车 正 代 表 二 者 之 间 的 质量 流量 .在 
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Fi F2.(12.2) Hh, RA SE u, v 便 可 以 看 作 是 物理 的 分 速 , 因为 只 差 
一 个 常数 因子 ， 这 时 RE 2， 于 是 方程 (13.2) 所 定义 的 流 画 数 
4 在 二 维 定常 无 旋 流 里 的 方程 便 是 
_ Ng wo NOY 

(1 n oe 2 oe +(1 “ss 0. (13.3) 

比较 方程 (13.3) 和 (9.11), 即 可 看 到 势 画 数 AU o WE 
同一 形式 的 方程 。 当 然 几 的 方程 和 由 的 方程 实际 的 内 容 是 不 相同 
的 ; u, o 和 中 之 间 的 关系 以 及 w,v 和 峭 之 间 的 关系 是 大 不 相同 
的 ， 例 如 方程 (13.1) 可 以 写成 


_r—l |u|? UD p 
(1 = isi) oe (13.4) 
这 里 面 的 四 是 对 应 于 驻 点 条 件 的 声速 ,而 方程 (12.8) 是 
PV 2 (OLY 4 (2y 
(Jz (24) + oy (13.5) 


自 方程 (9.4) 得 对 应 于 方程 (12.7) 的 关系 为 


at = (a? — I Jul’ (13.6) 


Fh F2( 9.11) F0(13.3 AEJK MERA DEB, BMRA, OA 
效 的 办 法 是 先 把 方程 线 化 ， 然 后 求解 线 化 方程 ， 把 方程 线 化 有 三 
种 方法 可 用 .第 一 种 方法 是 瑞 利 和 篇 森 〈Janzen) 提出 的 ,把 解 展 
为 Ma AAR, Mo 是 远离 物体 的 均匀 流 马 赫 数 ， BRR 
不 可 压 流 的 解 ,这 种 方法 在 Mo 很 小 时 是 适用 的 . 

第 二 种 方法 把 解 展 为 “厚度 比 ” 以 及 物体 轴线 与 远 痪 物体 均匀 
流 之 间 的 倾斜 角 的 升 需 级 数 。 厚 度 比 是 指 物体 在 主流 横断 方向 的 
厚度 与 物体 沿 主 流 方向 的 最 大 尺寸 之 比 。 这 种 方法 适用 于 亚 声 速 
或 超声 速 流 中 的 细 长 体 和 蒲 物体 . 其 一 级 近似 为 韶 劳 温 (Glaucrt)、 
普 朗 伦 和 阿 克 里 特 〈Ackeret) 等 所 发 展 ( 详 情 请 参看 第 VI 第 VU 
). 

第 三 种 线 化 的 方法 是 以 分 速 x 和 >” 或 以 等 价 的 六 二 |u| fu 
对 >z 轴 的 和 铬 角 8 作 自 变数 来 变换 方程 ， 这 就 是 速度 面 法 .以 VV 和 
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0 作 自 变数 的 势 画 数 由 的 方程 以 及 流 画 数 由 的 方程 可 以 直接 推导 
goo; 


图 A.13 


拿 相 邻 的 两 条 等 势 线 由 和 中 十 5 以 及 相 邻 的 两 条 流 线 少 和 
b+ op 所 组 成 的 网 格 来 看 ， 沿 等 势 线 距离 记 为 n, WREE 
为 ;。 设 流 线 的 曲率 和 牛 径 为 RC 图 A.13), 则 无 旋 条 件 是 


-E o, 
ðn R 
但 1/R = — 99/6s; & 
1 OV oe 
— — — — = 0 13. 
V On Os “ ( 7) 
连续 条 件 要 求 沙 上 4 与 48 十 ob 的 微 元 流 管 ， pV6n = const., Èp 
196 , 1 V 1 3n = 0 (13 8) 
p Os V Os ôn Os ° 
但 动力 学 方程 要 求 
pV OV 一 一 Op. 一 一 FRIA 
ðs ðs ðs 
所 以 
130, 1G yy i & 


p Os g Os V as? 
此 处 M 是 当地 马赫 数 V /a, 而 4 是 当地 声速 。 但 图 A. 13 又 表明 


1 06n _ 00 
n Os On” 
所 以 方程 (13.8) 变 成 


ôv , 30 
1—-m tX p E =o 13.9 
《 dF a ap e (13.9) 


现在 de = Vds, db = (0/0) V dn; 所 以 方程 (13.7) 和 (13.9) 
可 以 写成 


fil wv WwW o 13.10 
oY Op Op C ) 
及 
-mL 4 2 8 4g 13.11 
(1 > Oe + Op . C ) 
INV FD 6 VATE > A b HER ES, WF MIB KARR: 
ðv _ 1 ð$ 88 _ _ 1 Od 
Ob J 06° Ob J OV (13.12) 
V_ 1 db 88 _1 Ob 
Op J 00° ¢ J a’ 


此 处 了 是 牙科 上 毕 (Jacobian) 式 : 
j= alp., L) _ Ob Of _ Ob OL 


av, 9) OV 30 06 arv’ 


YF 13.12 MEA EE (13.10) 和 (13.11), 得 基本 的 速度 面 方程 
如 下 : 


2 1 Ob _ Ob 
oy æ av’? 
(13.13) 
26 (1 — m) LL Ob 
OV eV 00 


这 两 个 方程 现在 是 线性 的 了 ,因为 密度 p 和 马赫 数 M 都 只 是 了 了 的 
He. 从 这 两 个 方程 里 消去 上 或 $, MOR HAA, 这 
些 详 细 的 推导 见 第 VIE F Ga. 

轴 对 称 流 . 无 旋 轴 对 称 流 的 流 画 数 4 的 方程 可 以 从 方程 
(12.13),《12.16) 和 (12.17) 推 得 。 于 是 
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aw 


Op -2 ru, Od P ro, (13.14) 
Or o Ox p 
PY = (2EY 4 (2Y 13.15 
(5 ry lal = g Or (13.15) 
及 
(: -4)S4— uv _ Op ( £) oe 
a/ Ox? a OxOr a’/ Or? 
1 Ob o, (13.16) 
r Or 
RES SHARE 


2 2 v 2 v? 2 a 
Ce 2g A D 
(13.17) 
方程 (13.16) 和 (13.17) 又 不 是 线性 的 。 瑞 利 -篇 森 的 线 化 法 以 及 把 
解 展 为 厚度 参数 的 级 数 的 方法 在 这 里 仍 可 用 .不 好 的 是 ， 速 度 面 
变换 法 现在 不 能 把 $ 和 由 的 方程 线 化 了 . 


AL. 对 完全 气体 律 的 偏离 


在 前 几 节 的 讨论 中 常常 假定 气体 是 完全 气体 ， 然 而 有 许多 方 
程 是 写成 在 一 般 情况 下 也 成 立 的 形式 的 . 事实 上 ,只 要 方程 里 不 包 
括 cos cy 或 7 的话 , 那 就 是 一 般 的 方程 如果 包括 这 些 参 数 , 那 就 
是 适用 于 完全 气体 的 方程 ， 这 里 应 该 强调 指出 ， 流 体 可 以 是 理想 
的 ; 朗 无 粘性 又 不 传 热 , 但 却 不 是 完全 气体 , 完全 气体 是 必须 符合 
状态 方程 (5.7) 的 。 例 如 ,在 激 波 和 边界 层 之 外 , 粘性 和 热传导 的 
效应 都 可 以 略 去 不 计 , 然 而 假如 压强 十 分 高 ,温度 很 高 或 者 反之 温 
度 非常 低 , 那 末 气 体 的 性 质 是 和 完全 气体 的 性 质 有 明显 的 偏离 的 ， 
高 超声 风 洞 的 喷 管 中 的 流 劲 ,就 是 这 种 流动 的 一 个 例子 ,由 于 膨胀 
比 十 分 巨大 ， 喷 管 中 发 生 了 很 极端 的 情况 .这 种 偏离 和 由 有 一 定 
变化 率 的 过 程 所 产生 的 效应 完全 不 同 ， 郎 与 调整 到 热力 学 平衡 的 
过 程 所 产生 的 效应 完全 不 同 。 另 一 方面 ， 粘 性 和 热传导 就 是 以 有 
限 变 化 率 调整 到 热力 学 平衡 的 作用 结果 ， 因 而 朗 使 在 完全 气体 里 
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也 能 出 现 . 

有 两 种 偏离 完全 气体 律 的 情况 气体 可 以 很 准确 地 符合 完全 
气体 的 状态 方程 ,但 比 热 可 以 不 是 常数 ,这 种 气体 在 热学 上 是 完全 
的 ， 但 是 在 热量 关系 上 是 非 完全 的 气体 ， 或 者 气体 可 以 具有 常 比 
热 ,但 却 不 符合 完全 气体 的 状态 方程 ;这 和 神气 体 在 热量 关系 上 是 完 
全 的 ， 但 在 热学 上 是 非 完全 的 。 把 这 两 种 不 完全 加 以 区 分 是 很 有 
理 册 的。 如果 气 体 的 密度 不 大 ,分 子 之 闻 的 平均 距离 比较 大 ,在 这 
样 的 平均 距离 上 ,分 子 之 间 的 相互 作用 可 以 略 去 不 计 , 那 末 状 态 方 
程 就 是 完全 气体 的 状态 方程 。 另 一 方面 ， 等 容 比 热 是 为 分 子 的 各 
种 量子 状态 之 间 的 内 能 统计 分 布 所 决定 的 ,温度 高 于 或 低 于 室温 ， 
比 热 都 可 能 偏离 常数 值 相当 大 .在 十 分 低温 下 , 比 热 下 降 , 因 为 这 
时 分 子 的 旋转 自由 度 未 受 激 发 ;在 十 分 高 漫 下 , 比 热 上 升 , 因为 分 
子 的 振动 自由 度 被 激发 ， 所 以 气体 在 热量 关系 上 的 非 完 全 性 是 气 
体 分 子 内 部 特殊 结构 的 结果 ， 在 近 于 室温 的 温度 之 下 ， 没 有 这 些 
效应 ;但 如 果 气 体 的 密度 很 大 ,因而 分 子 之 间 的 平均 距 高 变 得 很 小 
时 ， 在 这 样 的 平均 距离 上 分 子 之 间 会 发 生 相 互 作 用 。 这 时 状态 方 
程 就 不 是 完全 气体 的 状态 方程 ， 气 体 的 热学 非 完全 性 是 由 于 分 子 
之 间 的 相互 作用 的 结果 ， 

可 以 借助 考虑 压强 为 零 的 气体 来 分 开 热量 关系 上 的 非 完 全 性 
和 热学 的 非 完全 性 ， 设 压强 为 零 时 的 内 能 为 co; AA eo 便 只 是 温 
RT A, RAD SRE 


eT, va) =|" [7 (22) e] dow + elt), Gaa 
式 中 wp 是 气体 的 比 容 , 部 


4, (14.2) 
p 


Usp 


括号 的 下 标 表示 该 参数 保持 不 变 ， —~HERM AS (ATR 
KA) 时 的 内 能 eT) 和 状态 方程 , 就 可 以 用 方程 (14.1 ) 计 算 在 任 
何 压强 ( 比 容 ) 及 温度 下 的 内 能 。 烩 只 是 
h = e + pos, (14.3) 
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对 于 完全 气体 来 说 ,方程 (14.1) 中 的 积分 为 需 , XN ARE RS 
压强 无 关 . 
等 类 过 程 无 热 加 入 .所 以 
0 = de + pdvsps 
AJEA), SES LG HE 


Op 
aT 工 (Sr). 
(2) 一 一 Ne, (14.4) 
dvsy/s dey 
dT 


TERS RERETREN SMR RY, SSR Ze, 
就 可 以 计算 伯 努 利 方程 (8.2) 中 的 积分 了 . 
声速 * 可 以 计算 如 下 : 


(的 -的 (后 
do/s P N dysp/s OT/ rsp \ dvsp/s 
+(3e) I. 
Ovsp T 


用 方程 (14.4) 代 (dT /dvsp),, 声速 平方 的 最 后 公式 是 


=v, OT /na 一 人 êL) (14.5) 
P deo Ovsp T ° 
dT 


如 果 是 完全 气体 ,方程 (14.5) 就 变 成 了 方程 (1.4) 那 样 的 简单 形式 ， 
声速 只 是 温度 的 一 个 男 数 . 但 方程 (14.5) 却 说 明 声 速 一 般 地 瞩 决 
定 于 气体 的 温度 又 决定 于 密度 的 . 

在 分 析 流 场 时 ,方程 (14.1) 和 (14.5) 可 以 提供 计算 占 重 要 地 位 
的 那些 物理 量 的 资料 ， 在 原则 上 用 这 两 个 方程 可 以 把 气体 偏离 完 
全 气体 律 的 影响 确定 出 来 。 实际 的 计算 过 程 当 然 是 很 繁 的 ， 得 用 
数值 法 和 图 解法 .例如 布 焉 曼 (Busemann)([201, 第 421 R) SHE 
出 过 怎样 构 作 计算 二 维 超声 速 流 所 需要 的 特征 线 图 。 因为 通常 偏 
离 完全 气体 并 不 大 ， 更 合 逻 辑 的 办 法 是 采用 以 完全 气体 的 解 作 为 
未 扰 解 的 摄 动 法 .这 种 办 法 ， 钱 学 森 在 研究 范 德 巨 耳 斯 气体 的 一 
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ETI Hs AIL ; 艾格 斯 (Eggers )P7) 处 理 同样 的 问题 也 用 过 ， 他 
用 的 是 贝 瑟 劳 特 〈Berthelot) 状态 方程 . 


A.15. 一 般 正 交 和 坐标 柔 中 的 运动 方程 的 表达 式 
REXER. 设 一 般 正 交 坐 标 系 上 的 一 点 坐标 为 w。 p> 
7， 在 该 点 上 在 a, 8, 7 增加 的 方向 上 三 微 元 线段 为 hida, hdl, 
Ady, F#EMK ds 是 
(ds)? = hilda) + RABY + RBar). 
如 果 01, 025 03 ERK V E a, 8, y 增 大 的 方向 上 的 分 量 , # 
VX V 的 三 分 量 为 二 n, 6, W 


v-Vv=—! | 之 
Ayhzh3 


2 Chhavi) H 一 一 sp Chaba) 


+ 一 一 人 | 
Fy Cavs) 


及 
一 awe Zo 303) 一 > Chwa) |, 
_ "一 二 [c Chwi) 一 2 - Chawa) ， 
一 po | 2 Chw) — 0 FCs). 
梯度 的 三 分 量 是 


18 1ð3 1ð 


Ay da” hy pb hs ðr’ 


应 力 张 量 的 各 分 量 可 以 写成 : 

Nae = — p + peua + (x — ta) (eaa + epp + err), 
xpe = —p + pepp + (w ras ) (ena + epg + err), 
xry = — p + uerr + (x 一 二 Cesa 十 epp + err), 


mpr 一 wyp = Lepr, Mra = ay 一 Mera, Wap 一 Npa 一 [Caps 


. 58 。 


式 中 各 ¢ Ala, b, y MARE u, v, 来 表示 , 则 为 


I ea = be 4 2 Ok X 8h 
2 Ay Oa Ayhz op hzhı Or 
1 ene = 2 4 w Oh ad Oho 
2 hy og hhz Oy hiha da” 


= 1 3w 4 4 3h y Z Dj 


= y = 一 - 一 一 


2 hor h Oa hh OB? 
em = 2 (2) + A OL (<), 
hy op hg ha Or 
er. =% (+) 4% 之 )， 
hs Or Ay hy ða ha 


eam 0 (2) Aa (4) 
hy Oa hy ha op Ay 


动力 学 方程 (2.7) 可 以 写成 
av LN EVD 
+(V. VV- a +y(+ 1v) 


Vx(VxXV)=xX+lv-n, 


APH X EAD ME, x EN DK, V- nr 的 三 个 分 量 是 : 
0 向 : 


1 [2 (hahsraa) + 一 一 sp Cheat) + 一 一 2. Caierroj 
Ayhohs 
pd hy, 1h 1 Oh 
hiha op hihs Or hih2 Oa 
_ an 1 3h, 
hihs Oc 
B: 
1 
TAn Loe Orte) + Chatna) + -2 Cerar) 
hyhoh3 
+ apy 2 22 op -1 Oh, 一 arr i hz 
hhz Or hh Oa hha og 
rl On 
hh OB 
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7 向 : 


1 [2 0 0 | 
一 一 一 | 一 一 (hhsrra) + —— (hsh 十 一 一 (Caporyy) 
hhh L Ox ¢ 2 ara) 65 《 3 pr) ay 《 14277YY 
Fe O01 
zhi Oa hah2 og hah1 Or 
~ rpg | Oh 
h3hz Or” 


FREH. 用 坐标 Ts Ws Xs 使 
xı = r cos w, x = rsinw, = x, 


如 取 r, w, x TBA a, B, To Ml A, = 1, h =r, h = 1, 于 是 


7 Ow Ox? 
g =L 0% Ov pa On _ Ov 
r Ow Ox” Ox Or ” 7 
=i 3 — 1 0 
ta r 8 (rez) r Oo ， 
(V .了 V)V 的 三 个 分 量 是 
r 向 : oe 十 了 au or 
Or 7 Ow Ox r 
co 向 : ze p 2 Ôe war we, . 
Or r Ow Ox r r’ 
xja: u Ow 4 2 2w w Oe 
Or r Ow Ox 
拉 普 拉 斯 算 子 V 是 
y= Z 19,1 0 © 
ðr? r ðr r? Ow? Ox? f 
定义 速度 势 为 
zx 一 过, 1 OH = 26 。 
Or r ðw Ox 
速度 势 的 方程 可 以 用 方程 (9.7) 写 出 来 . 7 
应 变 牵 的 各 分 量 是 . 
sa 60 >» 


1 = Ou 1 — 1 Ov “u _ Ow. 
Err D9? T fw = — O + —, Crr T 一; 
2 Or 2 r Ow > 2 Ox 
eo = 2 3e y Ov en = ow 
r ðw 0 Ox Or 
9 ( Z) 1 Ou 
ero = r 一 一 | 一 一 
Or r r Ow 
V zx 的 三 分 量 是 : 
r 向 : Orr 十 1 Oyo 十 On, 十 Wyp 一 Tow, 
ðr r Ow Ox r 


o fij: Ot y L Onow | Dror | 2 Hrs, 
Or r Ow Ox r 


z 向 : -ze 4 L ros | Ones | tar 
Or r Ow Ox r 
球 极 坐标 . 用 坐标 >, 0, w, 使 
xı = r sin cos w, x= rsinOsinw, x3 = r cos 0 
如 分 别 合 +, 0,0 Ha, p, y, W 
h =1, h=r, kh =rsinð, 
于 是 


VV =i Ê (rn) + -iĝ sino) +- -二 22 
E= = [2 Casing) 一 2e), 
EEE L 

E= COn) — a 


(V. V)V 的 三 分 量 是 : 
7 向 : OH p Or | w Du tu 


Or r 0 rsin Ow r ? 
0 向: Ov pLa w Ov | uv_ wv cod, 
Or r 00 rsin@ Ow r r? 


of: PCL A v Ow w Ow | wu + vw cot 0 


Or r 8 r sinb Ow r r 
拉 普 拉 斯 算 子 V? 是 
1 ð 0 1 0 。 0 
vw? = st LY ( 2 2) + _- ( 8 -二 
? or \" Or/) sind 00 00 
1 g 
r? sin? 0 ðw” 
速度 势 由 的 定义 是 
y= 06,106 1 de, 
Or r 06 rsin@ Ow 
速度 势 的 方程 可 以 根据 方程 (9.7) 求 出 . 
应 变 达 的 各 分 量 是 : 
1 1, Ll ee 
2 0 2 r 06 r 
i, =t Bw p” y veo, 
2 rsin@ Ow r r 


C gw 


r 00 rsin@ Oo” 


eu, = —L— ov + 2 (4) 
rsin Ow Or \r 


eg =r- (2) 4 1 


sin @ 


Or \r r 06° 
Y .的 三 分 量 是 
Br 1 Orr 1 On 1 
r : 一 二 十 一 一 一 十 一 于 十 一 (22, — ey 
向 Or r 00 rsin Ow r (2x ™ 


一 tgo + mp cot); 
0 向 : On, 十 工 aroe 4+—1 GETM 
Or r 06 rsinO Ow 
一 xww) cot 9 + 32915 


1 
十 -一 
r [Crog 


Orp 1 Omen 1 On 1 
op: Ste 4 + Cree 4 1 Ory igg 
Or r 00 rsin@ Ow r 3m 


+ 229. cot @), 
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定常 轴 对 称 流 的 到 性 坐标 柔 . 2: 代表 沿 子午 流 线 的 长 度 ， 

n 代表 在 子午 平面 内 垂直 于 流 线 的 长 度 ; 7 是 高 对 称 轴线 的 距离 . 
在 子午 面 内 的 速度 矢 V, 为 两 个 量 所 规定 : 一 个 是 子午 面 内 的 分 
RE om BDE Vn 相对 于 对 称 轴 的 余 角 9。 与 子午 面 垂直 
的 分 速度 是 w。 Es, n 和 圆周 方向 的 w 分 别 当 作 a, 8 和 7Y， 大 
注意 ,由 于 对 称 的 绿 故 , 只 有 对 ;和 7 的 导数 ， 于 是 
V-V = Pene. + Pusing + pon 2; 


> 
S 


g = 1 Ores) n= — 1 Gres) | t= _ Ovm + „30. 
r On r Os On Os 
应 变 率 的 各 分 量 是 : 
1 Ov,, 1 00 1 — Um 0 
es 一 > + Can S Um D> T Gow — ~~ sin, 
2 Os 2 ðn 2 r 
_,9 (+) _,9 (¢:) _ Ovum a0 
Caw = FTN j Cas = rl), Ens = + tm. 
On \r Os \r On Os 
动力 学 方程 是 
2 
PUm Oem _ p= sin @ = Or: 十 Orn + rs (4 sine + 2) 
Os r Os On r On 
+ aq (+ cose 一 2 æ) — Kan 8 
Os On 
一 ron > sin 6 + pX,; 
r 
2 
ov, 22. — po = cos0 = Orn + Or nn + x,, (4 sino +2 æ) 
Os 7 Os On r n 


+ tan + cose 一 tow cos + 0X,3 
r r 


pom O(re,) = Or + On pw + Kurs (2+ sin @ + æ) 
r ðs Os On r 


n 


+ aye (2+ cos0 — 2) + OX。。 
r 0 


定常 二 维 流 的 票 性 坐标 柔 . 全 :代表 沿 流 线 的 长 度 , > 代表 
与 之 垂直 的 长 度 。 速度 矢 V 为 其 数值 > RV 对 x 轴 的 倾角 6 所 
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规定 。 于 是 
V+ pV = + py OF, 
Os On 
t= Ov +o 30 , 
n Os 
应 变 率 的 各 分 量 是 
1 _ Ov 1 — OO 
Ess Zs — Cnn v 3 
2 Os 2 On 
Cns ðv +v 30 
On Os 
动力 学 方程 是 
y 22. w= .are 4 08 _ Qn, OF — wp, 28 + pX,, 
Os Os On On Os On 
po? 28. — Onon 十 Oran + 2x, .60 + pX 
On Os On On 


第 卡尔 坐标 系 . ARREARS MRD x, y, z， 相 应 


的 分 速 为 u, v, w。， 则 连续 方程 是 
o 4 Hou) , Ape) , lew) o 


Or Ox Oy Oz 
涡 舌 的 三 个 分 量 是 
gpa Ow _ Ov n= 04 _ Ow ca Oe Ou 
Oy Oz” Oz Ox” Ox Oy° 
粘性 应 力 各 分 量 是 
Ou G ) ( av 8w) 
Tee 一 24 — + 一 -十 一 -十 - 芝 )， 
Ox 四 fj\ar By Oz 


MIX, Y, Z 为 彻 体力 的 三 分 量 , 则 动力 学 方程 是 


Ou Ou Ou Ou Op 
Îi oy OE pv OH + pw OH = — L + ox 
fos or "Oy ”Bz Ox 


十 OT rx 十 ÔT ry 十 OTe 
Ox Oy Oz . 


Ov Ov Ov Ov _ oppy 


+ pu 2% + pv E + pw 
ae or OP ay ”Bz ay 


十 OT yx + OT yy + OT ys 
Ox Oy Oz 


Ow Ow Ow Ow Op 
= + pu —= + pu = + pw = — È 4 oZ 
Pa or Oy Bz ðs 


4 Tz 4 Oy: 十 charg 
Ox Oy Oz 


现在 耗 散 西数 @ 是 
ERORE) 


+ (Oe + Be) + (2 4 Bey) 
Oz Ox Ox Oy 


能 量 方程 是 
D _{ 9 0 0 2 ) 
“~~ (h =(2- 4,2 2 -2 
p © G “Br Oy ”Be h 


e 5 ¢ 
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